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Abstract. Now that the low temperature properties of quantum-mech- 
anical many-body systems (bosons) at low density, p, can be exam- 
ined experimentally it is appropriate to revisit some of the formulas 
deduced by many authors 4-5 decades ago, and to explore new regimes 
not treated before. For systems with repulsive (i.e. positive) interaction 
potentials the experimental low temperature state and the ground state 
are effectively synonymous - and this fact is used in all modeling. In 
such cases, the leading term in the energy /particle is 2-Kh 2 ap/m where 
a is the scattering length of the two-body potential. Owing to the del- 
icate and peculiar nature of bosonic correlations (such as the strange 
iV 7 / 5 law for charged bosons), four decades of research failed to estab- 
lish this plausible formula rigorously. The only previous lower bound 
for the energy was found by Dyson in 1957, but it was 14 times too 
small. The correct asymptotic formula has been obtained by us and 
this work will be presented. The reason behind the mathematical dif- 
ficulties will be emphasized. A different formula, postulated as late as 
1971 by Schick, holds in two dimensions and this, too, will be shown 
to be correct. With the aid of the methodology developed to prove the 
lower bound for the homogeneous gas, several other problems have been 
successfully addressed. One is the proof by us that the Gross-Pitaevskii 
equation correctly describes the ground state in the 'traps' actually used 
in the experiments. For this system it is also possible to prove complete 
Bose condensation and superfluidity, as we have shown. On the fron- 
tier of experimental developments is the possibility that a dilute gas in 
an elongated trap will behave like a one-dimensional system; we have 
proved this mathematically. Another topic is a proof that Foldy's 1961 
theory of a high density Bose gas of charged particles correctly describes 
its ground state energy; using this we can also prove the N 7 ^ 5 formula 
for the ground state energy of the two-component charged Bose gas pro- 
posed by Dyson in 1967. All of this is quite recent work and it is hoped 
that the mathematical methodology might be useful, ultimately, to solve 
more complex problems connected with these interesting systems. 
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Foreword 

At the conference "Perspectives in Analysis" at the KTH, Stockholm, 
June 23, 2003, one of us (E.H.L.) contributed a talk with the title above. 
This talk covered material by all the authors listed above. This contribution 
is a much expanded version of the talk and of |L4j . It is based on, but 
supersedes, the article [Lj3SY . 
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1. Introduction 

Schrodinger's equation of 1926 defined a new mechanics whose Hamilton- 
ian is based on classical mechanics, but whose consequences are sometimes 
non-intuitive from the classical point of view. One of the most extreme cases 
is the behavior of the ground (= lowest energy) state of a many-body system 
of particles. Since the ground state function ^/{xi, ...,a?jv) is automatically 
symmetric in the coordinates Xj G K 3 of the ./V particles, we are dealing 
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necessarily with 'bosons '. If we imposed the Pauli exclusion principle (anti- 
symmetry) instead, appropriate for electrons, the outcome would look much 
more natural and, oddly, more classical. Indeed, the Pauli principle is es- 
sential for understanding the stability of the ordinary matter that surrounds 
us. 

Recent experiments have confirmed some of the bizarre properties of 
bosons close to their ground state, but the theoretical ideas go back to 
the 1940's - 1960's. The first sophisticated analysis of a gas or liquid of 
interacting bosons is due to Bogolubov in 1947. His approximate theory as 
amplified by others, is supposed to be exact in certain limiting cases, and 
some of those cases have now been verified rigorously (for the ground state 
energy) — 3 or 4 decades after they were proposed. 

The discussion will center around five main topics. 

1. The dilute, homogeneous Bose gas with repulsive interaction (2D 
and 3D). 

2. Repulsive bosons in a trap (as used in recent experiments) and the 
"Gross-Pitaevskii" equation. 

3. Bose condensation and superfluidity for dilute trapped gases. 

4. One-dimensional behavior of three-dimensional gases in elongated 
traps. 

5. Foldy's "jellium" model of charged particles in a neutralizing back- 
ground and the extension to the two-component gas. 

Note that for potentials that tend to zero at infinity 'repulsive' and 'pos- 
itive' are synonymous — in the quantum mechanical literature at least. In 
classical mechanics, in contrast, a potential that is positive but not mono- 
tonically decreasing is not called repulsive. 

The discussion below of topic 1 is based on |LYlj and |LY2j . and of topic 2 
on jLSeYlj and |LSeY2j . See also [LY3l ILSe Y3I ISe2l ILSe Y4j . 

The original references for topic 3 are |LSej and |LSeY5j . but for trans- 
parency we also include here a section on the special case when the trap is 
a rectangular box. This case already contains the salient points, but avoids 
several complications due the inhomogeneity of the gas in a general trap. 
An essential technical tool for topic 3 is a generalized Poincare inequality, 
which is discussed in a separate section. 

Topic 4 is a summary of the work in LS eY6j . 

The discussion of topic 5 is based on |LSoj and |LSo2 . 

Topic 1 (3-dimensions) was the starting point and contains essential ideas. 
It is explained here in some detail and is taken, with minor modifications 
(and corrections), from |LY3j . In terms of technical complexity, however, 
the fifth topic is the most involved and can not be treated here in full detail. 

The interaction potential between pairs of particles in the Jellium model 
in topic 5 is the repulsive, long-range Coulomb potential, while in topics 
1-4 it is assumed to be repulsive and short range. For alkali atoms in the 
recent experiments on Bose Einstein condensation the interaction potential 
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has a repulsive hard core, but also a quite deep attractive contribution of 
van der Waals type and there are many two body bound states |PSj . The 
Bose condensate seen in the experiments is thus not the true ground state 
(which would be a solid) but a metastable state. Nevertheless, it is usual to 
model this metastable state as the ground state of a system with a repulsive 
two body potential having the same scattering length as the true potential, 
and this is what we shall do. In this paper all interaction potentials will be 
positive. 

2. The Dilute Bose Gas in 3D 

We consider the Hamiltonian for N bosons of mass m enclosed in a cubic 
box A of side length L and interacting by a spherically symmetric pair 
potential v(\xi — Xj\): 

N 

R N = -/^ A *+ Yl v(\xi-Xj\). (2.1) 

i=l l<i<j<N 

Here X{ £R 3 , i = 1,...,N are the positions of the particles, Aj the Lapla- 
cian with respect to Xj, and we have denoted h 2 /2m by p for short. (By 
choosing suitable units p could, of course, be eliminated, but we want to 
keep track of the dependence of the energy on Planck's constant and the 
mass.) The Hamiltonian (j2.1j) operates on symmetric wave functions in 
L 2 (A N , dx\ ■ ■ ■ dxjv) as is appropriate for bosons. The interaction poten- 
tial will be assumed to be nonnegative and to decrease faster than 1/r 3 at 
infinity. 

We are interested in the ground state energy Eq(N,L) of (|2.1j) in the 
thermodynamic limit when N and L tend to infinity with the density p = 
iV/L 3 fixed. The energy per particle in this limit is 

e (p)= lim E (pL\L)/(pL 3 ). (2.2) 

Our results about eo(p) are based on estimates on Eq(N, L) for finite N and 
L, which are important, e.g., for the considerations of inhomogeneous sys- 
tems in |LSeYlj . To define Eq(N, L) precisely one must specify the boundary 
conditions. These should not matter for the thermodynamic limit. To be 
on the safe side we use Neumann boundary conditions for the lower bound, 
and Dirichlet boundary conditions for the upper bound since these lead, 
respectively, to the lowest and the highest energies. 

For experiments with dilute gases the low density asymptotics of eo(p) is 
of importance. Low density means here that the mean interparticle distance, 
p-i/ 3 i s much larger than the scattering length a of the potential, which is 
defined as follows. The zero energy scattering Schrodinger equation 



-2pAip + v(r)rp = 



(2.3) 
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has a solution of the form, asymptotically as | x \ = r — > oo (or for all r > Rq 
if v(r) = for r > Ro), 

Mx) = l~a/\x\ (2.4) 
(The factor 2 in ()2.3|) comes from the reduced mass of the two particle 
problem.) Writing il)o{x) = uo(|^|)/|x| this is the same as 

a = Hm r - ^M, (2.5) 
r-*oo Uo{f) 

where uq solves the zero energy (radial) scattering equation, 

-2/iug(r) + v(r)u {r) = (2.6) 

with uq(0) = 0. 

An important special case is the hard core potential v(r) = oo if r < a 
and v(r) = otherwise. Then the scattering length a and the radius a are 
the same. 

Our main result is a rigorous proof of the formula 

e Q (p) w Airfipa (2.7) 

for pa 3 <C 1, more precisely of 

Theorem 2.1 (Low density limit of the ground state energy). 

hm ^Bi = 1. (2.8) 
pa 3 ^o Air /i pa 

This formula is independent of the boundary conditions used for the defi- 
nition of eo(p) • It holds for every positive radially symmetric pair potential 
such that v(r)r 2 dr < oo for some R, which guarantees a finite scattering 
length, cf. Appendix A in |LY2j . 

The genesis of an understanding of e$(p) was the pioneering work |Boj of 
Bogolubov, and in the 50's and early 60's several derivations of ()2.8|) were 
presented jHY], EE], 

even including higher order terms: 
^ = 1 + ^{ P a^ + 8 (f - VS) (pa") log( P a 3 ) + 0(pa 3 ) (2.9) 



These early developments are reviewed in |L2j . They all rely on some special 
assumptions about the ground state that have never been proved, or on the 
selection of special terms from a perturbation series which likely diverges. 
The only rigorous estimates of this period were established by Dyson, who 
derived the following bounds in 1957 for a gas of hard spheres |Dlj : 



1 e (p) l + 2Y 1 /3 

< ow < Z (2.10) 



WV2 ~ 4nppa ~ (1 - yi/3)2 

with Y = 4npa 3 /3. While the upper bound has the asymptotically correct 
form, the lower bound is off the mark by a factor of about 1/14. But for 
about 40 years this was the best lower bound available! 

Under the assumption that ()2.8|) is a correct asymptotic formula for the 
energy, we see at once that understanding it physically, much less proving it, 
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is not a simple matter. Initially, the problem presents us with two lengths, 
a <C p -1 / 3 at low density. However, ((2.8(1 presents us with another length 
generated by the solution to the problem. This length is the de Broglie wave- 
length, or 'uncertainty principle' length (sometimes called 'healing length') 

4~(H~ 1/2 - (2-11) 

The reason for saying that l c is the de Broglie wavelength is that in the 
hard core case all the energy is kinetic (the hard core just imposes a ip = 
boundary condition whenever the distance between two particles is less than 
a). By the uncertainty principle, the kinetic energy is proportional to an 
inverse length squared, namely i c . We then have the relation (since pc? is 
small) 

a < p- 1/3 < 4 (2.12) 

which implies, physically, that it is impossible to localize the particles rela- 
tive to each other (even though p is small). Bosons in their ground state are 
therefore 'smeared out' over distances large compared to the mean particle 
distance and their individuality is entirely lost. They cannot be localized 
with respect to each other without changing the kinetic energy enormously. 

Fermions, on the other hand, prefer to sit in 'private rooms', i.e., l c is 
never bigger than p" 1 / 3 by a fixed factor. In this respect the quantum 
nature of bosons is much more pronounced than for fermions. 

Since ((2.8|) is a basic result about the Bose gas it is clearly important 
to derive it rigorously and in reasonable generality, in particular for more 
general cases than hard spheres. The question immediately arises for which 
interaction potentials one may expect it to be true. A notable fact is that 
it is not true for all v with a > 0, since there are two body potentials with 
positive scattering length that allow many body bound states. (There are 
even such potentials without two body bound states but with three body 
bound states |Ba| .) For such potentials (|2.8|) is clearly false. Our proof, 
presented in the sequel, works for nonnegative v, but we conjecture that 
1(2.8(1 holds if a > and v has no TV-body bound states for any N. The 
lower bound is, of course, the hardest part, but the upper bound is not 
altogether trivial either. 

Before we start with the estimates a simple computation and some heuris- 
tics may be helpful to make ((2.8(1 plausible and motivate the formal proofs. 

With ipo the zero energy scattering solution, partial integration, using 
(IP) and (23J), gives, for R > R , 

/ {2/x|V'0o| 2 + v\^o\ 2 }dx = 8irpa ( 1 — — ) — > Snpa for R — > oo. 
J\x\<R v R J 

(2.13) 

Moreover, for positive interaction potentials the scattering solution mini- 
mizes the quadratic form in ((2.13() for each R > Rq with the boundary con- 
dition V'oG^I = R) = (1 — a/R). Hence the energy Eq(2, L) of two particles 
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in a large box, i.e., L 3> a, is approximately 87rpa/L 3 . If the gas is suffi- 
ciently dilute it is not unreasonable to expect that the energy is essentially 
a sum of all such two particle contributions. Since there are N(N — l)/2 
pairs, we are thus lead to E (N, L) « A-npaN(N — 1)/L 3 , which gives (|2.8|) 
in the thermodynamic limit. 

This simple heuristics is far from a rigorous proof, however, especially for 
the lower bound. In fact, it is rather remarkable that the same asymptotic 
formula holds both for 'soft' interaction potentials, where perturbation the- 
ory can be expected to be a good approximation, and potentials like hard 
spheres where this is not so. In the former case the ground state is approxi- 
mately the constant function and the energy is mostly potential: According 
to perturbation theory Eq(N, L) ~ N(N — 1)/(2L 3 ) J v(\x\)dx. In particular 
it is independent of p, i.e. of Planck's constant and mass. Since, however, 
J v(\x\)dx is the first Born approximation to 8irp,a (note that a depends on 
p\), this is not in conflict with (|2.8[) . For 'hard' potentials on the other hand, 
the ground state is highly correlated, i.e., it is far from being a product of 
single particle states. The energy is here mostly kinetic, because the wave 
function is very small where the potential is large. These two quite different 
regimes, the potential energy dominated one and the kinetic energy dom- 
inated one, cannot be distinguished by the low density asymptotics of the 
energy. Whether they behave differently with respect to other phenomena, 
e.g., Bose-Einstein condensation, is not known at present. 

Bogolubov's analysis |Boj presupposes the existence of Bose-Einstein con- 
densation. Nevertheless, it is correct (for the energy) for the one-dimensional 
delta- function Bose gas |LL| . despite the fact that there is (presumably) no 
condensation in that case |PiStj . It turns out that BE condensation is not 
really needed in order to understand the energy. As we shall see, 'global' 
condensation can be replaced by a 'local' condensation on boxes whose size 
is independent of L. It is this crucial understanding that enables us to prove 
Theorem 1.1 without having to decide about BE condensation. 

An important idea of Dyson was to transform the hard sphere potential 
into a soft potential at the cost of sacrificing the kinetic energy, i.e., effec- 
tively to move from one regime to the other. We shall make use of this 
idea in our proof of the lower bound below. But first we discuss the simpler 
upper bound, which relies on other ideas from Dyson's beautiful paper |Dlj . 

2.1. Upper Bound. The following generalization of Dyson's upper bound 
holds |LSeYlj . |SeT] : 

Theorem 2.2 (Upper bound). Let pi = (N-l)/L 3 andb= (A-Kpi/2,)- 1 / 3 . 
For non-negative potentials v and b > a the ground state energy of \2.1\) with 
periodic boundary conditions satisfies 

E (N,L)/N < 4vr m a b ^ Kbl 8 2 ■ (2.14) 

I 1 ~~ f ) 
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Thus in the thermodynamic limit (and for all boundary conditions) 

-W < (2.15) 



infipa ~ (1 - Y 1 / 3 ) 
provided Y = Anpo? /3 < 1. 

Remark. The bound (|2.14|) holds for potentials with infinite range, provided 
b > a. For potentials of finite range Ro it can be improved for b > Rq to 

1 _ f^ 2 + I f«) 3 
£ (iV, L)/N < Air^a ^ -U^-. (2.16) 

(!-!) 

Proof. We first remark that the expectation value of (|2,lj) with any trial wave 
function gives an upper bound to the bosonic ground state energy, even if 
the trial function is not symmetric under permutations of the variables. The 
reason is that an absolute ground state of the elliptic differential operator 
(|2.1j) (i.e., a ground state without symmetry requirement) is a nonnegative 
function which can be symmetrized without changing the energy because 
(|2.1j) is symmetric under permutations. In other words, the absolute ground 
state energy is the same as the bosonic ground state energy. 
Following |Dlj we choose a trial function of the following form 

*(x x , . . . ,x N ) = F 1 (x 1 ) ■ F 2 (x 1 ,x 2 ) ■ --Fn^i, . . .,x N ). (2.17) 

More specifically, F\ = 1 and F{ depends only on the distance of Xi to its 
nearest neighbor among the points x*x,..., Xj„i (taking the periodic bound- 
ary into account): 

Fi(xx, . . . ,Xi) = f(ti), U = mm(\xi - Xj\, j = 1, . . . ,i - 1) , (2.18) 
with a function / satisfying 

< / < 1, /' > 0. (2.19) 

The intuition behind the ansatz (|2.17|) is that the particles are inserted into 
the system one at the time, taking into account the particles previously in- 
serted. While such a wave function cannot reproduce all correlations present 
in the true ground state, it turns out to capture the leading term in the en- 
ergy for dilute gases. The form (|2.17|) is computationally easier to handle 
than an ansatz of the type Ili<j/(l^» — ^j'Dj which might appear more 
natural in view of the heuristic remarks after Eq. (|2.13|) . 
The function / is chosen to be 

/(r) = //o(r)//o(6) for0<r<6, 

JK ' [1 forr>6, V 7 

with fo(r) = uo(r)/r the zero energy scattering solution defined by Q2.6|) . 
The estimates (|2.14(l and (j2.16|) are obtained by somewhat lengthy compu- 
tations similar as in |Dlj . but making use of (|2.13|) . For details we refer to 
jLSeYlj and jSeTj. ■ 
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2.2. Lower Bound. It was explained previously in this section why the 
lower bound for the bosonic ground state energy of (|2.1j) is not easy to 



obtain. The three different length scales (|2.12|) for bosons will play a role in 
the proof below. 

• The scattering length a. 

• The mean particle distance p -1 / 3 . 

• The 'uncertainty principle length' £ c , defined by pi~ 2 = eo(p), i.e., 
4 ~ (pa)- 1 / 2 . 

Our lower bound for eo(p) is as follows. 

Theorem 2.3 (Lower bound in the thermodynamic limit). For a 

positive potential v with finite range and Y small enough 

i°M_>(l_(7Yi/i7) (2.21) 
Anppa 

with C a constant. If v does not have finite range, but decreases faster than 
1/r 3 (more precisely, v(r)r 2 dr < oo for some R) then an analogous 
bound to 12.21]) holds, but with CY 1 / 17 replaced by o(l) as Y — ► 0. 

It should be noted right away that the error term —CY 1 / 17 in (J2.2H) is 
of no fundamental significance and is not believed to reflect the true state 
of affairs. Presumably, it does not even have the right sign. We mention in 
passing that C can be taken to be 8.9 jSelj . 

As mentioned at the beginning of this section after Eq. (|2.2[) . a lower 
bound on Eq(N,L) for finite iV and L is of importance for applications to 
inhomogeneous gases, and in fact we derive (|2.21|) from such a bound. We 
state it in the following way: 

Theorem 2.4 (Lower bound in a finite box). For a positive potential v 
with finite range there is a 5 > such that the ground state energy of \2. 1\) 
with Neumann boundary conditions satisfies 



E (N, L) /N > Airppa (l-C Y 1/17 ) (2.22 



for all N and L with Y < 5 and L/a > C'Y~ 6 / 17 . Here C and C are 
positive constants, independent of N and L. (Note that the condition on L/a 
requires in particular that N must be large enough, N > (const. )Y~ 1 t 17 .) 
As in Theorem \2J-ft such a bound, but possibly with a different error term 
holds also for potentials v of infinite range that decrease sufficiently fast at 
infinity. 

The first step in the proof of Theorem 12.41 is a generalization of a lemma 
of Dyson, which allows us to replace v by a 'soft' potential, at the cost of 
sacrificing kinetic energy and increasing the effective range. 

Lemma 2.5. Let v(r) > with finite range Rq. Let U(r) > be any 
function satisfying f U(r)r 2 dr < 1 and U(r) =0 for r < Rq. Let 8cl 3 
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be star shaped with respect to (e.g. convex with £ B). Then for all 
differentiable functions ip 

[ M |VVf + i*#| 2 ] >^a [ U\^\ 2 . (2.23) 

JB 

Proof. Actually, (|2.23|) holds with fi\ Vip(x)\ 2 replaced by the (smaller) radial 
kinetic energy, fi\dip(x) / dr\ 2 , and it suffices to prove the analog of (|2.23f) for 
the integral along each radial line with fixed angular variables. Along such 
a line we write ip(x) = u(r)/r with u(0) = 0. We consider first the special 
case when U is a delta- function at some radius R> Rq, i.e., 

U(r) = ±5(r-R). (2.24) 

For such U the analog of (|2.23j) along the radial line is 

{fj,[u'(r) - (u(r)/r)] 2 + \v{r)\u{r)] 2 }dr > j ° , m|2/D2 'f. ^ ^ „ 

I fia\u{R)\ j R it it < R\ 

(2.25) 

where R\ is the length of the radial line segment in B. The case R\ < R is 
trivial, because [i\dij) / dr\ 2 + ^IV 7 ! 2 — 0- (Note that positivity of v is used 
here.) If R < R\ we consider the integral on the left side of (|2,25|) from 
to R instead of R\ and minimize it under the boundary condition that 
u(0) = and u(R) is a fixed constant. Since everything is homogeneous in 
u we may normalize this value to u(R) = R — a. This minimization problem 
leads to the zero energy scattering equation (|2.6[) . Since v is positive, the 
solution is a true minimum and not just a stationary point. 

Because v(r) =0 for r > Rq the solution, uq, satisfies uo(r) = r — a for 
r > Rq. By partial integration, 

R 

{^[u' Q {r)-{uQ{r)/r)] 2 + \v{r)\uQ{r)] 2 }dr = fia\R-a\/R > fia\R-a\ 2 /R 2 . 
o 

(2.26) 

But \R — a\ 2 /R 2 is precisely the right side of (|2.25|) if u satisfies the normal- 
ization condition. 

This derivation of (|2.23|) for the special case (|2.24|) implies the general 
case, because every U can be written as a superposition of ^-functions, 
U(r) = j R~ 2 5{r-R) U(R)R 2 dR, and / U{R)R 2 dR < 1 by assumption. ■ 

By dividing A for given points x\, . . . , xjy into Voronoi cells Bi that contain 
all points closer to X{ than to Xj with j ^ i (these cells are star shaped w.r.t. 
Xi , indeed convex) , the following corollary of Lemma 12.51 can be derived in 
the same way as the corresponding Eq. (28) in 1)1. 

Corollary 2.6. For any U as in Lemma \2. 5\ 

H N > fiaW (2.27) 
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with W the multiplication operator 

N 

W(x 1 ,...,x N )=^U(t l ), (2.28) 
i=i 

where ti is the distance of to its nearest neighbor among the other points 
xj, j = 1,...,N, i.e., 

U(xi, . . . , xjy) = min \xi — xA. (2.29) 

(Note that ti has here a slightly different meaning than in (|2.18jl . where it 
denoted the distance to the nearest neighbor among the Xj with j < i — 1.) 

Dyson considers in |Dlj a one parameter family of C/'s that is essentially 
the same as the following choice, which is convenient for the present purpose: 

UR{r) = (3(« 3 " HI)" 1 *» *><r<R (2 . 30) 
I otherwise. 

We denote the corresponding interaction (|2.28j) by Wr. For the hard core 
gas one obtains 

E(N, L) > sup inf fiaW R (xi, . . . , x N ) (2.31) 

R (X!,...,X N ) 

where the infimum is over (xi,... ,xn) G A n with \x{ — Xj\ > Ro = a, 
because of the hard core. At fixed R simple geometry gives 

inf, W R (£ 1} ...,x N )>(^-^) (2.32) 
{xi,...,s N ) \R 6 pR b J 

with certain constants A and B. An evaluation of these constants gives 
Dyson's bound 

E(N,L)/N>^A^pa. (2.33) 

The main reason this method does not give a better bound is that R 
must be chosen quite big, namely of the order of the mean particle distance 
p -1 / 3 , in order to guarantee that the spheres of radius R around the iV 
points overlap. Otherwise the infimum of Wr will be zero. But large R 
means that Wr is small. It should also be noted that this method does not 
work for potentials other than hard spheres: If \x^ — xA is allowed to be less 
than Ro, then the right side of (|2.31j) is zero because U(r) = for r < Rq. 

For these reasons we take another route. We still use Lemma 12.51 to get 
into the soft potential regime, but we do not sacrifice all the kinetic energy 
as in (|2.27jl . Instead we write, for e > 

H N = eH N + (1 - e)H N > eT N + (1 - e)H N (2.34) 

with Tjy = — Yli ^-i an d use ()2.27|) only for the part (1 — e)Hn- This gives 

H N > eT N + (1 - e)fiaW R . (2.35) 
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We consider the operator on the right side from the viewpoint of first order 
perturbation theory, with eTjv as the unperturbed part, denoted Hq. 

The ground state of Hq in a box of side length L is *f?o(xi, • • • >%n) = 
L~ 3N / 2 and we denote expectation values in this state by (-)o- A computa- 
tion, cf. Eq. (21) in [LYT] (see also Eqs. {EI5H32D)), gives 

47rp(l-i) > (W R ) /N 

> 4,p (1 - (1 - M) 3 (1 + 4vrp(i? 3 - i? 3 )/3) 1 . 

(2.36) 

The rationale behind the various factors is as follows: (1 — ■h) comes from 
the fact that the number of pairs is N(N - 1) /2 and not N 2 /2, (1 - 2R/L) 3 
takes into account the fact that the particles do not interact beyond the 
boundary of A, and the last factor measures the probability to find another 
particle within the interaction range of the potential U R for a given particle. 

The estimates (|2.36|) on the first order term look at first sight quite promis- 
ing, for if we let L —> 00, N — » 00 with p = N/L 3 fixed, and subsequently 
take R — ► 0, then (W R )o/N converges to 4irp, which is just what is desired. 
But the first order result (|2.36|) is not a rigorous bound on Eq(N,L), we 
need error estimates, and these will depend on e, R and L. 

We now recall Temple 's inequality [T] for the expectation values of an op- 
erator H = Hq + V in the ground state (-)o of Hq. It is a simple consequence 
of the operator inequality 

(H — Eq){H — E\) > (2.37) 
for the two lowest eigenvalues, Eq < E±, of H and reads 

Eo > (H) ~ (2-38) 

provided E\ — {H)q > 0. Furthermore, if V > we may use E l > E{ 0) = 

second lowest eigenvalue of Hq and replace E\ in (|2.38|) by Ey . 
From 1)2.36(1 and (|2.38J) we get the estimate 

E (N,L) 



N 



>4irpap(l-£(p,L,R,e)) (2.39) 



with 



l-£(p,L,R,e) 

= (1-^(1-^) (l-¥) 3 (l + fP^ 3 -^))" 
x L _ ^((W 2 R ) - (W R )l) \ 
X \ (W r ) {e{ 0) - na(W R )o) J ' 
To evaluate this further one may use the estimates (|2.36|) and the bound 



(2.40) 



(W 2 R ) < 3 R3 * R3 (W R ) (2.41) 
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which follows from = 3(i? 3 — Rq)~ 1 Ur together with the Schwarz in- 
equality. A glance at the form of the error term reveals, however, that it 
is not possible here to take the thermodynamic limit L — > oo with p fixed: 
We have E^ = eir^/i/L 2 (this is the kinetic energy of a single particle in 
the first excited state in the box), and the factor e[ ^ — pa(Wn)o in the de- 
nominator in 1)2. 40j) is, up to unimportant constants and lower order terms, 
~ (eL~ 2 — ap 2 L 3 ). Hence the denominator eventually becomes negative and 
Temple's inequality looses its validity if L is large enough. 

As a way out of this dilemma we divide the big box A into cubic cells of 
side length £ that is kept fixed as L — ► oo. The number of cells, L 3 /£ 3 , on the 
other hand, increases with L. The N particles are distributed among these 
cells, and we use (j2.4Ujl . with L replaced by £, N by the particle number, n, 
in a cell and p by n/£ 3 , to estimate the energy in each cell with Neumann 
conditions on the boundary. For each distribution of the particles we add 
the contributions from the cells, neglecting interactions across boundaries. 
Since v > by assumption, this can only lower the energy. Finally, we 
minimize over all possible choices of the particle numbers for the various 
cells adding up to N. The energy obtained in this way is a lower bound to 
Eq(N,L), because we are effectively allowing discontinuous test functions 
for the quadratic form given by H^. 

In mathematical terms, the cell method leads to 

E (N, L)/N > (pi 3 )- 1 inf Y, CnE (n, I) (2.42) 

n>0 

where the infimum is over all choices of coefficients c n > (relative number 
of cells containing exactly n particles), satisfying the constraints 

J>„ = 1, Y,c n n = p£ 3 . (2.43) 

n>0 n>0 

The minimization problem for the distributions of the particles among 
the cells would be easy if we knew that the ground state energy Eq(u,£) (or 
a good lower bound to it) were convex in n. Then we could immediately 
conclude that it is best to have the particles as evenly distributed among 
the boxes as possible, i.e., c n would be zero except for the n equal to the 
integer closest to p£ 3 '. This would give 

Eo{ ^ L) > Anpap (1 - £(p, £, R, e)) (2.44) 

i.e., replacement of L in (|2.39l) by £, which is independent of L. The blow 
up of £ for L — > oo would thus be avoided. 

Since convexity of Eo(n,£) is not known (except in the thermodynamic 
limit) we must resort to other means to show that n = 0(p£ 3 ) in all boxes. 
The rescue comes from superadditivity of Eq(u,£), i.e., the property 

E (n + n',£)>E (n,e)+E (n',£) (2.45) 
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which follows immediately from v > by dropping the interactions between 
the n particles and the n' particles. The bound (|2.45|) implies in particular 
that for any n,p£N with n> p 

77 

E (n, £) > [n/p] E (p, £) > —E (p, £) (2.46) 

2p 

since the largest integer [n/p] smaller than n/p is in any case > n/(2p). 

The way ()2.46j) is used is as follows: Replacing L by £, N by n and p by 
n/£ 3 in ()2.39[) we have for fixed R and e 

Ea(n,£) > ^p n (n-l)K(n,l) (2.47) 

with a certain function K(n,£) determined by (|2.4(H) . We shall see that K 
is monotonously decreasing in n, so that if p S N and n < p then 

E Q (n,£) > ^p n (n-l)K(p,£). (2.48) 

We now split the sum in (|2.42l) into two parts. For n < p we use (|2.48(l . and 
for n > p we use (|2.46|) together with (j2,48j) for n = p. The task is thus to 
minimize 

c n n(n - 1) + \ c n n(p - 1) (2.49) 

n<p n>p 

subject to the constraints (|2.43|) . Putting 

k := p£ 3 and t := ^ c n n < k (2.50) 

we have Xm>p c « n = k — t, and since n(n — 1) is convex in n and vanishes 
for 77 = 0, and Z~2n<p c n — 1) ^ ne expression 1)2. 49() is 

>t(t-l) + l(k-t)(p-l). (2.51) 

We have to minimize this for 1 < t < k. If p > 4£i the minimum is taken at 
t = k and is equal to k(k — 1). Altogether we have thus shown that 

E o(N, L) ^ ^ p f 1 _J_\ K{Ap e A) , (2 .5 2) 



N ~ ^ r \ p£ 3 / 

What remains is to take a closer look at K(4p£ 3 ,£), which depends on 
the parameters e and R besides £, and choose the parameters in an optimal 
way. From (|2.40j) and (|2.41j) we obtain 

KM = (l-s)(l-^) 3 (l + f(R s -Rl)y 1 

( 3 an \ 

X V 1 ~ k (R 3 ~ R 3 )(ire£~ 2 - ia£~Mn -I)))' 

The estimate (|2.47|) with this K is valid as long as the denominator in the 
last factor in (|2.53|) is > 0, and in order to have a formula for all n we can 
take as a trivial lower bound in other cases or when 1)2. 47|) is negative. 
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As required for (|2.48jl . K is monotonously decreasing in n. We now insert 
n = Ap£ 3 and obtain 

K(4p£ 3 ,£) > (l- e )(l-2K) 3 (l + (const.)Y{£/a) 3 {R 3 - R^/e 3 )' 1 

1 _ _P (const.)Y \ 

(R? - R 3 ) (e(a/£) 2 - (const.)Y 2 (£/a) 3 ) J 1 ' 

with Y = 4irpa 3 /3 as before. Also, the factor 

1 - \ ] = (1 - (const. )Y- 1 {a/£) 3 ) (2.55) 

in (|2.52|) (which is the ratio between nin — 1) and n 2 ) must not be forgotten. 
We now make the ansatz 

e~Y a , a/£~Y?, {R 3 - R 3 )/£ 3 ~ Y"< (2.56) 

with exponents a, j3 and 7 that we choose in an optimal way. The conditions 
to be met are as follows: 

• e(a/£) 2 — (const. )Y 2 {£/a) 3 > 0. This holds for all small enough Y, 
provided a + 5/3 < 2 which follows from the conditions below. 

• a > in order that e — > for Y — > 0. 

• 3/3 - 1 > in order that Y- 1 (a/£) 3 -> for for Y -> 0. 

• 1-3/3 + 7 > in order that Y (£ / a) 3 (R 3 - R 3 ) / £ 3 -» for for y 0. 

• 1— a — 2/3 — 7>0to control the last factor in (|2.54j) . 
Taking 

a = 1/17, (3 = 6/17, 7 = 3/17 (2.57) 

all these conditions are satisfied, and 

a = 3/3-l = l-3/? + 7 = l- a-2/3-7 = 1/17. (2.58) 

It is also clear that 2R/£ ~ y^/ 3 = y 1 / 17 , up to higher order terms. This 
completes the proof of Theorems 12.31 and 12.41 for the case of potentials 
with finite range. By optimizing the proportionality constants in (j2.56|) one 
can show that C = 8.9 is possible in Theorem 12.31 SelJ. The extension to 
potentials of infinite range but finite scattering length is obtained by approx- 
imation by finite range potentials, controlling the change of the scattering 
length as the cut-off is removed. See Appendix A in |LY2j and Appendix 
B in |LSeYlj for details. We remark that a slower decrease of the potential 
than 1/r 3 implies infinite scattering length. ■ 

The exponents 1)2. 5 7|) mean in particular that 

a < R « p- 1 ' 3 < £ < {pa)- 1 ' 2 , (2.59) 

whereas Dyson's method required R ~ p~ 1 ' 3 as already explained. The 
condition p^ 1 ' 3 <C £ is required in order to have many particles in each box 
and thus n(n — 1) ~ n 2 . The condition £ <C (pa)- 1 ' 2 is necessary for a 
spectral gap 3> eo(p) in Temple's inequality. It is also clear that this choice 
of £ would lead to a far too big energy and no bound for eo(p) if we had 
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chosen Dirichlet instead of Neumann boundary conditions for the cells. But 
with the latter the method works! 

3. The Dilute Bose Gas in 2D 

In contrast to the three-dimensional theory, the two-dimensional Bose gas 
began to receive attention only relatively late. The first derivation of the 
correct asymptotic formula was, to our knowledge, done by Schick L S^ for a 
gas of hard discs. He found 

e(p)^47r W \ln(pa 2 )\-\ (3.1) 

This was accomplished by an infinite summation of 'perturbation series' 
diagrams. Subsequently, a corrected modification of [Sj was given in [HFM . 
Positive temperature extensions were given in |Poj and in |FH| . All this work 
involved an analysis in momentum space, with the exception of a method 
due to one of us that works directly in configuration space |Llj . Ovchinnikov 
[01 derived (jHHJ) by using, basically, the method in |Llj . These derivations 
require several unproven assumptions and are not rigorous. 

In two dimensions the scattering length a is defined using the zero energy 
scattering equation 1)2.3)1 but instead of ip{r) ~ 1 — a/r we now impose the 
asymptotic condition ip(r) ~ ln(r/a). This is explained in the appendix to 

ra. 

Note that in two dimensions the ground state energy could not possibly 
be eo(p) ~ 47r/x/9o as in three dimensions because that would be dimen- 
sionally wrong. Since eo(p) should essentially be proportional to p, there is 
apparently no room for an a dependence — which is ridiculous! It turns out 
that this dependence comes about in the ln(pa 2 ) factor. 

One of the intriguing facts about (|3.1|) is that the energy for N parti- 
cles is not equal to iV(iV — l)/2 times the energy for two particles in the 
low density limit — as is the case in three dimensions. The latter quan- 
tity, Eq(2,L), is, asymptotically for large L, equal to 87rpL~ 2 [ln(L 2 /a 2 )] 
(This is seen in an analogous way as (|2.13|) . The three-dimensional bound- 
ary condition -00 (1^1 = R) = 1 — a/R is replaced by V'od^l = R) = ln(R/a) 
and moreover it has to be taken into account that with this normalization 
\\ipo\\ 2 = (volume)(ln(i2/a)) 2 (to leading order), instead of just the volume in 
the three-dimensional case.) Thus, if the N(N—l)/2 rule were to apply, 1)3.1)1 
would have to be replaced by the much smaller quantity Airpp [ln(L 2 /a 2 )] 
In other words, L, which tends to oo in the thermodynamic limit, has to 
be replaced by the mean particle separation, p -1 / 2 in the logarithmic fac- 
tor. Various poetic formulations of this curious fact have been given, but 
the fact remains that the non-linearity is something that does not occur in 
more than two dimensions and its precise nature is hardly obvious, physi- 
cally. This anomaly is the main reason that the two-dimensional case is not 
a trivial extension of the three-dimensional one. 
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Eq. (|3.1j) was proved in |LY2j for nonnegative, finite range two-body po- 
tentials by finding upper and lower bounds of the correct form, using similar 
ideas as in the previous section for the three-dimensional case. We discuss 
below the modifications that have to be made in the present two-dimensional 
case. The restriction to finite range can be relaxed as in three dimensions, 
but the restriction to nonnegative v cannot be removed in the current state 
of our methodology. The upper bounds will have relative remainder terms 
0(| ln^a 2 )!^ 1 ) while the lower bound will have remainder 0(| \n(pa 2 )\~ 1 ^). 
It is claimed in |HFM| that the relative error for a hard core gas is negative 
and 0(ln| ln(/oa. 2 ) 1 1 ln(/3a 2 )| _1 ), which is consistent with our bounds. 

The upper bound is derived in complete analogy with the three dimen- 
sional case. The function /o in the variational ansatz (|2.2()j) is in two dimen- 
sions also the zero energy scattering solution — but for 2D, of course. The 
result is 

Eo{N > L)/N * HbM-*P* (1 + 0([ln(6/a)r 1} ) • (3 - 2) 

The minimum over b of the leading term is obtained for b = (27rp) -1 / 2 . 
Inserting this in (j3,2[) we thus obtain 

E (N, L)/N < (1 + 0(| ln^a 2 )!" 1 )) . (3.3) 

To prove the lower bound the essential new step is to modify Dyson's 
lemma for 2D. The 2D version of Lemma 12.51 is: 

Lemma 3.1. Let v(r) > and v(r) = for r > Rq. Let U(r) > be any 
function satisfying 

/•oo 

/ U(r)\n(r/a)rdr < 1 and U(r) = for r < R . (3.4) 
Jo 

Let B C M. 2 be star-shaped with respect to (e.g. convex with G B). Then, 
for all functions rp in the Sobolev space H 1 (B), 

{ijl\V^{x)\ 2 + \v(\x\)\^(x)\ 2 ) dx>fif U{\x\)\ip{x)\ 2 dx. (3.5) 



Proof. In polar coordinates, r, 6, one has (V^l 2 > \dil)/dr\ 2 . Therefore, it 
suffices to prove that for each angle 9 € [0, 2ir), and with ip(r,9) denoted 
simply by /(r), 

R(fi) rR{9) 

(fi\df{r)/dr\ 2 + \v{r)\f{r)\ 2 )rdr> fi U(r)\f(r)\ 2 rdr, (3.6) 

Jo 

where R{0) denotes the distance of the origin to the boundary of B along 
the ray 6. 

If R{0) < -Ro then (|3.6|) is trivial because the right side is zero while the 
left side is evidently nonnegative. (Here, v > is used.) 

If R(0) > Rq for some given value of 0, consider the disc T>(9) = {x G 
M 2 : < \x\ < R(9)} centered at the origin in M? and of radius R(9). Our 
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function / defines a spherically symmetric function, x i— ► /(|x|) on T>{8), 
and (|3.6j) is equivalent to 

HV/(|f|)| 2 + ^(l^l)l/(l^l)| 2 l dx>J U(\x\)\f{\m 2 dx. 



V(0) \ 2 / JV{0) 

(3.7) 

Now choose some R E {Ro, R{0)) and note that the left side of (|3.7|) is 
not smaller than the same quantity with T>{9) replaced by the smaller disc 
T^R = {x £ M 2 : < \x\ < R}. (Again, v > is used.) We now minimize 
this integral over T>r, fixing f{R). This minimization problem leads to the 
zero energy scattering equation. Plugging in the solution and integrating by 
parts leads to 



2vr J U\df{r)/dr\' + -v{r)\f{r)\ 2 j rdr > ^-^\f{R)\\ (3.8) 

The proof is completed by multiplying both sides of (|3.8|) by U {R)R \n{R/a) 
and integrating with respect to R from Rq to R{9). ■ 

As in Corollary 12.61 Lemma 13.11 can be used to bound the many body 
Hamiltonian Hn from below, as follows: 

Corollary 3.2. For any U as in Lemma \3.1\ and any < e < 1 

H N > eT N + (1 - e)nW (3.9) 
with Tjv = —fi Y2i=i an d 

W{x\, . . . , xjf) = / U ( mm \xi — Xj\. ) . (3.10) 
i=i K ' 

For U we choose the following functions, parameterized by R > Rq: 

UR(r) = {-(«)-' for«,<r<fl 
I otherwise 

with ^(ii) chosen so that 

C/ R (r) ln(r/o)r dr = 1 (3.12) 

Ro 

for all R> Rq, i.e., 

= / m{r/a)r dr = \ {R 2 {m{R 2 /a 2 ) - l) - R 2 {ln{R 2 /a 2 ) - l) } . 
JRo 

(3.13) 

The nearest neighbor interaction (|3.1()|) corresponding to Ur will be denoted 
W R . 

As in Subsection 12.21 we shall need estimates on the expectation value, 
{Wr)q, of Wr in the ground state of eT/v of (|3.9j) with Neumann boundary 
conditions. This is just the average value of Wr in a hypercube in M. 2N . 
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Besides the normalization factor v(R), the computation involves the volume 
(area) of the support of Ur, which is 

A(R) = tt(R 2 - Rl), (3.14) 

In contrast to the three-dimensional situation the normalization factor 
v{R) is not just a constant {R independent) multiple of A{R); the fac- 
tor ln(r/a) in ()3.4|) accounts for the more complicated expressions in the 
two-dimensional case. Taking into account that Ur is proportional to the 
characteristic function of a disc of radius R with a hole of radius Ro, the 
following inequalities for n particles in a box of side length £ are obtained 
by the same geometric reasoning as lead to ()2.36j) . cf. |LYlj : 



(Wr) > 
(W R ) < 



u(R) 



\- 2 -fY i-(i-Q)(«-D 



n 



u{R) 



i-(i-Q) 



(n-l) 



(3.15) 
(3.16) 

Q = A(R)/£ 2 (3.17) 

being the relative volume occupied by the support of the potential Ur. Since 
Ur = v{R)~ 1 Ur we also have 



with 



R/0 



< 



n 



v(R) 



{Wr) . 



As in [EyT| we estimate [1 - (1 - Q)^ 1 )] by 



This gives 



(n-l)Q> [l - (1 - Q) (n_1) 
n(n — 1) 



> 



(n-l)Q 
l + (n-l)Q 



(Wr) > 



(Wr) < 



Q 



v(R) l + (n-l)Q' 
n(n — 1) 



v{R) 



Q 



From Temple's inequality fTj we obtain like in (|2.38|) the estimate 

K(wi)o - (w R )i) \ 



Eo(n,£) > (l-e)(W R ) 1 



(Wr) {e[ 0} -fi(W R ) ) 



where 



E 



(o) 



(3.18) 

(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



is the energy of the lowest excited state of eT n . This estimate is valid for 
/ n > (Wr)q, i.e., it is important that £ is not too big. 
Putting (|3.2U|) - (J3.22j) together we obtain the estimate 



, . n(n-l) A(R) 
E (n, I) > V ' —rJr K(n) 



£ 2 



v{R) 



(3.24) 
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with 

(1-25)2 / \ 

= <* - e > ■ i + (» - i)Q ' I 1 - ( , ^/g) - - 1) Q ) (3 ' 25> 

Note that Q depends on £ and R, and K depends on £, R and e besides n. 
We have here dropped the term (W/?)o in the numerator in (|3.22j) . which is 
appropriate for the purpose of a lower bound. 

We note that K is monotonically decreasing inn, so for a given n we may 
replace K(n) by K(p) provided p > n. As explained in the previous section, 
(|2.45|) - (|2.52jl . convexity of n *— ► n(n — 1) together with superadditivity of 
Eo(n, I) in n leads, for p = <ip£ 2 , to an estimate for the energy of iV particles 
in the large box when the side length L is an integer multiple of £: 

P A(R)( A \\ VIA _^ 



E (N, L)/N > {l-- 2 ) K(4p£') (3.26) 

with p = N/L 2 . 

Let us now look at the conditions on the parameters e, R and £ that have 
to be met in order to obtain a lower bound with the same leading term as 
the upper bound (|3.3|) . 

From (|3.13j) we have 

W) - - 1) (' " ^ r ' r °» < 3 ' 27 > 

We thus see that as long as a < R < p^ 1 ! 2 the logarithmic factor in the 
denominator in (|3.27|) has the right form for a lower bound. Moreover, for 
Temple's inequality the denominator in the third factor in Q3.25JI must be 
positive. With n = Ap£ 2 and v{R) > (const. )R 2 h\(R 2 /a 2 ) for R > R , this 
condition amounts to 

(const.)eln(i? 2 /a 2 )/^ 2 > p 2 f. (3.28) 
The relative error terms in ()3.26j) that have to be <C 1 are 

r 1 R nN 2 ^ f3 2Q^ 

£ ' p£ 2 ' V 9 ' eR 2 HR 2 /a 2 )- [ ' 

We now choose 

e~\]n(pa 2 T 1/5 , t ~ p'^lHpa 2 )^ 10 , R ~ ,cr 1/2 | ln(pa 2 )|~ 1/10 

(3.30) 

Condition (|3.28|) is satisfied since the left side is > (const.) | ln(pa 2 )| 3//5 and 
the right side is ~ | ln(pa 2 )| 2 / 5 . The first three error terms in (|3.29f) are all 
of the same order, | hi(pa 2 )\~ 1 ^ , the last is ~ | ln(/9a 2 )| _1 / 5 (ln | m(/9a 2 )|) _1 . 
With these choices, (j3.26j) thus leads to the following: 

Theorem 3.3 (Lower bound). For all N and L large enough such that 
L > (const. )p~ l,2 \ InOoa 2 )! 1 / 10 and N > (const.)] \n(pa 2 )\ 1/5 with p = N/L 2 , 
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the ground state energy with Neumann boundary condition satisfies 

E (N, L)/N > (l - 0(| Hpa 2 )\'^)) . (3.31) 

In combination with the upper bound (|3.3|) this also proves 
Theorem 3.4 (Energy at low density in the thermodynamic limit). 

lim 1 i?^ 2 ,| i =l (3-32) 

P a 2 ^o 4irpp\ ln^a 2 )!" 1 

where eo(p) = \\m.^—, O0 EQ{N,p^ 1 ^ 2 N 1 / 2 )/N. This holds irrespective of 
boundary conditions. 



As in the three-dimensional case, Theorem 13.41 is also valid for an infinite 

R 



range potential v provided that v > and for some R we have v(r)r dr < 



oo, which guarantees a finite scattering length. 

4. Generalized Poincare Inequalities 

This section contains some lemmas that are of independent mathematical 
interest, but whose significance for the physics of the Bose gas may not be 
obvious at this point. They will, however, turn out to be important tools 
for the discussion of Bose-Einstein condensation (BEC) and superfluidity in 
the next section. 

The classic Poincare inequality |LLoj bounds the L 9 -norm of a function, 
/, orthogonal to a given function g in a domain /C, in terms of some L p -norm 
of its gradient in 1C. For the proof of BEC we shall need a generalization 
of this inequality where the estimate is in terms of the gradient of / on a 
subset £1 C K, and a remainder that tends to zero with the volume of the 
complement Vl c = K, \ £l. For superfluidity it will be necessary to generalize 
this further by adding a vector potential to the gradient. This is the most 
complex of the lemmas because the other two can be derived directly from 
the classical Poincare inequality using Holder's inequality. The first lemma 
is the simplest variant and it is sufficient for the discussion of BEC in the 
case of a homogeneous gas. In this case the function g can be taken to 
be the constant function. The same holds for the second lemma, which 
will be used for the discussion of superfluidity in a homogeneous gas with 
periodic boundary conditions, but the modification of the gradient requires a 
more elaborate proof. The last lemma, that will be used for the discussion of 
BEC in the inhomogeneous case, is again a simple consequence of the classic 
Poincare and Holder inequalities. For a more comprehensive discussion of 
generalized Poincare inequalities with further generalizations we refer to 
PeY7| . 

Lemma 4.1 (Generalized Poincare inequality: Homogeneous case). 

Let K, C M 3 be a cube of side length L, and define the average of a function 
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/ G L\K) by 



(f)K = j3 J^f(x)dx . 



There exists a constant C such that for all measurable sets C K, and all 
f G i/ 1 (/C) the inequality 



[ \f{x)-{f)^\ 2 dx<cU 2 f \Vf(x)\ 2 dx + \n c \ 2 / 3 [ |V/( 



x)\ 2 dx 



holds. Here Vt c = fC\£l, and \ ■ \ denotes the measure of a set. . 

Proof. By scaling, it suffices to consider the case L = 1. Using the usual 
Poincare-Sobolev inequality on K. (see |LLoj . Thm. 8.12), we infer that there 
exists a C > such that 

\\f - {f)K\\hvo < hc\\vf\\ 2 L6/H!C) 

< ^(llV/|li 6/6(n) + ||V/|| 2 L6/5(nc) ) . (4.2) 
Applying Holder's inequality 

[|V/[| L 6/B (n) < IIV/H^^IOI 1 ^ 

(and the analogue with O replaced by Q c ), we see that (|4.1j) holds. ■ 

In the next lemma tC is again a cube of side length L, but we now replace 
the gradient V by 

V v := V + *(0,0, p/L), (4.3) 

where ip is a real parameter, and require periodic boundary conditions on 
JC. 

Lemma 4.2 (Generalized Poincare inequality with a vector poten- 
tial). For any \ip\ < tt there are constants c > and C < oo such that 
for all subsets Q, C K, and all functions f G H l (K,) with periodic boundary 
conditions on K, the following estimate holds: 

2 

||V<p/||£a (n j > -p H/Hi 2 (/C) + p 11/ - (/kllz^/c) 

-c(||V v /||i 2(K) + ^11/11^) (^) V2 • (4.4) 

i^ere |J7 C | is i/ie volume of Q c = K, \ O, i/ie complement of £1 in tC. 

Proof. We shall derive (|4.4|) from a special form of this inequality that holds 
for all functions that are orthogonal to the constant function. Namely, for 
any positive a < 2/3 and some constants c > and C < oo (depending only 
on a and \<p\ < tt) we claim that 



V^lli, { n)> ^11*00-5(^1 IIV^II 2 ^) , (l.o) 
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provided {l,h)/c = 0. (Remark: Eq. (|4.5|) holds also for a = 2/3, but the 
proof is slightly more complicated in that case. See |LSeY7j .) If (|4.5|) is 
known the derivation of 1)4. 4JI is easy: For any /, the function h = f — 
L~ 3 (l,f)jc is orthogonal to 1. Moreover, 

\\v<pHh<M = llv^llWi - liv^f 



v u HL 2 (n) — II w <p n \\L 2 (K) ~ II v v' t ilL2(n c ) 

h m -^\(L-^,f)K\ 2 (l + ^ 



and 



+2^Re(L-^ 2 ,f) lc (V v f,L-^ 2 ) nc 
< \\^f\\h m -^\(L~ 3/2 ,f)d 2 

+WwilWzll/«W>) 



^+ C „,„2 V 2 



1/2 



(4.6) 



L 2 



= ^(ll/lll^-K^/ki 2 ) 



+^H/-i:-3 (1)/k ||2 i2(jc) (47) 



Setting a = 1, using ||V^/t||£2nn < || ^ ipfWtfQC) i n the last term in 1)4.5)1 
and combining (|4~5|) . l{P|t and (|4~7|) gives (jOjl with C = \<p\ + C. 

We now turn to the proof of (|4.5|) . For simplicity we set L = 1. The 
general case follows by scaling. Assume that (j4.5|) is false. Then there 
exist sequences of constants C n — > oo, functions h n with ||/in||L 2 (A:) = 1 an d 
{l,h n )/c = 0, and domains O n C /C such that 

\\h{K)}<^- ( 4 - 8 ) 

We shall show that this leads to a contradiction. 

Since the sequence h n is bounded in L 2 (fC) it has a subsequence, denoted 
again by h n , that converges weakly to some h £ L 2 (fC) (i.e., (g,h n )/c — > 
(5, for all 5 £ L 2 (/C)). Moreover, by Holder's inequality the L P (S7^) norm 
IIV^/tnllipfQc) = (f Qlk \SJ lf ,h(x)\ p dx) 1 l' p is bounded by \&n\ a/2 \\V <phn\\lJ>{K.) 
for p = 2/(a + 1). From (j4.8)) we conclude that llV^/inHi^r^ ) is bounded 
and also that HV^/inlliP^) < ll^ 7 v>^n||L 2 (r2 n ) is bounded. Altogether, V ' v h n 
is bounded in L p (/C), and by passing to a further subsequence if necessary, 
we can therefore assume that V v /i n converges weakly in L p (/C). The same 
applies to Vh n . Since p = 2/(a + 1) with a < 2/3 the hypotheses of the 
Rellich-Kondrashov Theorem |LLo) Thm 8.9] are fulfilled and consequently 
h n converges strongly in L 2 {KL) to h (i.e., \\h — h n \\ L 2^ — > 0). We shall now 
show that 

liminf ||V^/i n ||| 2( n n ) > \N<ph\\h(K) • ( 4 - 9 ) 
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This will complete the proof because the h n are normalized and orthogonal 
to 1 and the same holds for h by strong convergence. Hence the right side 
of (|4.9j) is necessarily > ip 2 , since for \ip\ < n the lowest eigenvalue of — V£, 
with constant eigenfunction, is non-degenerate. This contradicts ()4.8|) . 

Eq. (|4.9|) is essentially a consequence of the weak lower semicontinuity 
of the L 2 norm, but the dependence on Q n leads to a slight complica- 
tion. First, Eq. (|4.8|) and C n — > oo clearly imply that — > 0, be- 
cause HV^/inll^^) > <p 2 - By choosing a subsequence we may assume that 

J2 n l^nl < 00 ■ For some fixed N let fi^r = K. \ U n >N^n- Then Qn C fi n for 
n > N. Since HV^/inU^^ \ is bounded, V v h n is also bounded in L 2 (f27v) 
and a subsequence of it converges weakly in L 2 (Qn) to V v h. Hence 

hminf ||V^/i„||| 2{nn) > liminf ||V v M£a ( n„) > \\V v h\\% { n N) • ( 4 - 10 ) 

Since fijv C ^at+i and Un^n = £ (up to a set of measure zero), we can 
now let N — > oo on the right side of (|4.1Uj) . By monotone convergence 
this converges to IIV^/iH^^)- This proves (|4.9f) which, as remarked above, 
contradicts (|4.8j) . ■ 

The last lemma is a simple generalization of Lemma 14.11 with /C C M. m 
a bounded and connected set that is sufficiently nice so that the Poincare- 
Sobolev inequality (see |LLol Thm. 8.12]) holds on fC. In particular, this is 
the case if fC satisfies the cone property |LLoj (e.g. if K, is a rectangular 
box or a cube). Moreover, the constant function on K, is here replaced by a 
more general bounded function. 

Lemma 4.3 (Generalized Poincare inequality: Inhomog. case). For 

d > 2 let K, C W 1 be as explained above, and let h be a bounded function 
with Jr h = 1 . There exists a constant C ( depending only on K, and h ) such 
that for all measurable sets Q C K, and all f G H l {K) with j^fhdx = 0, 
the inequality 

J^\f(x)\ 2 dx<C^\Vf(x)\ 2 dx+^ (4.11) 

holds. Here \ ■ \ denotes the measure of a set, and f2 c = fC \ fL 

Proof. By the usual Poincare-Sobolev inequality on K, (see |LLol 
Thm. 8.12]), 



< C||V/|| 

< 2C (||V/||2 2d/(d+2)(n) + \\Vf\\ 2 L2d/(d+2)m ) , (4.12) 
if d > 2 and J l<c fh = 0. Applying Holder's inequality 

\\Vf\\ L 2d /{ d + 2 ){n) < [|v/[| L3( n ) |n| 1/d 
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(and the analogue with Q replaced by O c ), we see that (j4.4jl holds with 
C = 2\)C\ 2 / d C. U 

5. Bose-Einstein Condensation and Superfluidity for 
Homogeneous Gases 

5.1. Bose-Einstein Condensation. Bose-Einstein condensation (BEC) is 
the phenomenon of a macroscopic occupation of a single one-particle quan- 
tum state, discovered by Einstein for thermal equilibrium states of an ideal 
Bose gas at sufficiently low temperatures [E|. We are here concerned with 
interacting Bose gases, where the question of the existence of BEC is highly 
nontrivial even for the ground state. Due to the interaction the many body 
ground state is not a product of one-particle states but the concept of a 
macroscopic occupation of a single state acquires a precise meaning through 
the one-particle density matrix. Given the normalized ground state wave 
function this is the operator on L 2 (M. d ) (d = 2 or 3) given by the kernel 

7(2, x) = N j ^ (x,X)^ (x,X)dX , (5.1) 

where we introduced the short hand notation 

N 

X = (a?2, . . . , xn) and dX = dxj. (5-2) 



Then J j(x, x)dx = Tr [7] = N. BEC in the ground state means, by defini- 
tion, that this operator has an eigenvalue of order N in the thermodynamic 
limit. Since 7 is a positive kernel and, hopefully, translation invariant in the 
thermodynamic limit, the eigenfunction belonging to the largest eigenvalue 
must be the constant function L~ d l 2 . Therefore, another way to say that 
there is BEC in the ground state is that 

J j(x,y)dxdy = 0(N) (5.3) 

as N — > 00, L — > 00 with N/L d fixed; more precisely Eq. (j5,3[) requires that 
there is a c > such that the left side is > cN for all large N. This is also 
referred to as off- diagonal long range order. Unfortunately, this is something 
that is frequently invoked but has so far never been proved for many body 
Hamiltonians with genuine interactions — except for one special case: hard 
core bosons on a lattice at half-filling (i.e., N = half the number of lattice 
sites) . The proof is in |KLSj and |DLSj . 

The problem remains open after more than 75 years since the first inves- 
tigations on the Bose gas [Bl [E] • Our construction in Section |2] shows that 
(in 3D) BEC exists on a length scale of order p _1 / 3 y -1 / 17 which, unfortu- 
nately, is not a 'thermodynamic' length like volume 1 / 3 . The same remark 
applies to the 2D case of Section 3, where BEC is proved over a length scale 
p-^lHpa 2 )] 1 / 10 . 
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In a certain limit, however, one can prove (|5 .,'-?[) . as has been shown in 
|LSej . In this limit the interaction potential v is varied with N so that the 
ratio a/L of the scattering length to the box length is of order 1/JV, i.e., the 
parameter Na/L is kept fixed. Changing a with N can be done by scaling, 
i.e., we write 

v(\x\) = ± Vl (\x\/a) (5.4) 

for some v\ having scattering length 1, and vary a while keeping v\ fixed. 
It is easily checked that the v so defined has scattering length a. It is 
important to note that, in the limit considered, a tends to zero (as iV~ 2 / 3 
since L = {N/p) 1 / 3 ~ N 1 / 3 for p fixed), and v becomes a hard potential 
of short range. This is the opposite of the usual mean field limit where the 
strength of the potential goes to zero while its range tends to infinity. 

We shall refer to this as the Gross-Pitaevskii (GP) limit since Na/L will 
turn out to be the natural interaction parameter for inhomogeneous Bose 
gases confined in traps, that are described by the Gross-Pitaevskii equation 
discussed in Sections and [7| Its significance for a homogeneous gas can 
also be seen by noting that Na/L is the ratio of pa to 1/L 2 , i.e., in the GP 
limit the interaction energy per particle is of the same order of magnitude as 
the energy gap in the box, so that the interaction is still clearly visible, even 
though a —* 0. Note that pa 3 ~ iV~ 2 in the GP limit, so letting N —* oo with 
p fixed and Na/L fixed can be regarded as a simultaneous thermodynamic 
and low density limit. For simplicity, we shall here treat only the 3D case. 

Theorem 5.1 (BEC in a dilute limit). Assume that, as N —* oo, p = 
N/L 3 and g = Na/L stay fixed, and impose either periodic or Neumann 
boundary conditions for H . Then 

N™ooN~~h 1 1 7(f ' y^ d y = 1 • ( 5 - 5 ) 

The reason we do not deal with Dirichlet boundary conditions at this 
point should be clear from the discussion preceding the theorem: There 
would be an additional contribution ~ 1 /L 2 to the energy, i.e. of the same 
order as the interaction energy, and the system would not be homogeneous 
any more. Dirichlet boundary conditions can, however, be treated with the 
methods of Section 

By scaling, the limit in Theorem 15. II is equivalent to considering a Bose 
gas in a fixed box of side length L = 1, and keeping Na fixed as N — » oo, i.e., 
a ~ 1/N. The ground state energy of the system is then, asymptotically, iVx 
4irNa, and Theorem l5.1l implies that the one-particle reduced density matrix 
7 of the ground state converges, after division by N, to the projection onto 
the constant function. An analogous result holds true for inhomogeneous 
systems as will be discussed in Section 

The proof of Theorem 15 . 1 1 has two main ingredients. One is localization of 
the energy that is stated as Lemma 15.21 below. This lemma is a refinement 
of the energy estimates of Section 2.2 and says essentially that the kinetic 
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energy of the ground state is concentrated in a subset of configuration space 
where at least one pair of particles is close together and whose volume tends 
to zero as a — > 0. The other is the generalized Poincare inequality, Lemma 
14. II from which one deduces that the one particle density matrix is approxi- 
mately constant if the kinetic energy is localized in a small set. 

The localization lemma will be proved in a slightly more general version 
that is necessary for Theorem 15.11 namely with the gradient V replaced by 
= V + i(0,0,(p/L), cf. Eq. I)4.3jl . We denote by H' N the corresponding 
many-body Hamiltonian 1)2.1)) with in place of V. This generalization 
will be used in the subsequent discussion of superfluidity, but a reader who 
wishes to focus on Theorem 15.11 only can simply ignore the <p and the ref- 
erence to the diamagnetic inequality in the proof. We denote the gradient 
with respect to x\ by Vi, and the corresponding modified operator by Vi i¥ ,. 

Lemma 5.2 (Localization of energy). Let K, be a box of side length L. 
For all symmetric, normalized wave functions \P(a?i, . . . ,x*n) with periodic 
boundary conditions on K,, and for N > Y~ l l 17 , 



±{%H' N V) > (1 - const. Y 1 / 17 ) 

x ( A-Kfipa + \x \ dX I <ixi|Vi (/ ,\I'(xi,X)| 



where X = (x*2, • • • , x^), dX = Y[f=2 d%j> o- n d 



(5.6) 



&x = : mm \x\ — xj\ > R > (5-7) 

with R = aY- 5 / 17 . 

Proof. Since is symmetric, the left side of (|5.6j) can be written as 
/ dX [ dx 1 L|Vi^^(xi,X)| 2 + iy?;(|f 1 -x i |)|^(x 1 ,X)| 2 l . (5.8) 

JK*- 1 JK 1 j>2 

For any e > and R > this is 

> eT + (1 - e) (T m + /) + (!- e)T° ut , (5.9) 



with 



T = nf dX [ dfi|Vi|^(fi,X)|| 2 , (5.10) 
T in = Li[ dX [ dfi|Vi|*(fi,X)|| 2 , (5.11) 



2 



T° ut =/W dX dxiVi,^(fi,X) , (5.12 
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1=1 [ dX f dxi 'yv(\xi-x j \)\^f(xi,X)\ 2 . (5.13) 
Jkk-i Jk jzi 



i>2 



ft c ^ = {xi : \x\ — Xj\ < R for some j > 2} 



(5.14) 

is the complement of $7^, and the diamagnetic inequality |LLoj \^7tpf(x)\ 2 > 
| V|/(x)| | 2 has been used. The proof is completed by using the estimates used 
for the proof of Theorem 2.4, in particular 1)2.521) and (|2.54|) - (|2.57|) . which 
tell us that for e = Y 1 / 17 and R = aY' 5 ^ 7 



sT+ (1 - e)(T in + I) > (l - const. Y 1/17 )4ir[ipa 
as long as N > Y" 1 / 17 . 



(5.15) 



Proof of Theorem \5.1\ We combine Lemma fS/il fwith ip = and hence = 
Hn) with Lemma 14.11 that gives a lower bound to the second term on the 
right side of (|5.6|) . We thus infer that, for any symmetric with (^f,^>) = 1 
and for N large enough, 

±{*,H N *)(1- const. Y 1 /!?)- 1 



> 47r/xpa - CY 



U 2 iV 



L Jk n - x Jk 



(5.16) 



where we used that \U C \ < ^-NR 3 



const. L 3 Y 2 / 17 . Since the kinetic 
energy, divided by N, is certainly bounded independent of N, as the upper 
bound (|2.14|) shows, and since the upper and the lower bound to Eq agree 
in the limit considered, the positive last term in (|5.16|) has to vanish in the 
limit. I.e., we get that for the ground state wave function of 



lim / dX I dxi ^ (xi,X) - L- 3 [Ldx^ (x,X)] 

N^oo JfcN-i Jfc v 

This proves (|5.5|) . since 

/ dX ( dxi VoixuX) - L- 3 [Ldx^ (x,X)] 
Jk n -^ Jk 

7(x, x')dxdx' . 







(5.17) 



NL 3 



(5.18) 



KxK 
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5.2. Superfluidity. The phenomenological two-fluid model of superfluidity 
(see, e.g., |TT| ) is based on the idea that the particle density p is composed 
of two parts, the density p s of the inviscid superfluid and the normal fluid 
density p n . If an external velocity field is imposed on the fluid (for instance 
by moving the walls of the container) only the viscous normal component 
responds to the velocity field, while the superfluid component stays at rest. 
In accord with these ideas the superfluid density in the ground state is often 
defined as follows |HoM| : Let Eq denote the ground state energy of the 
system in the rest frame and E' the ground state energy, measured in the 
moving frame, when a velocity field v is imposed. Then for small v 

f = ^ + (p s // 3 )imv 2 + 0(|v| 4 ) (5.19) 

where A" is the particle number and m the particle mass. At positive tem- 
peratures the ground state energy should be replaced by the free energy. 
(Remark: It is important here that (|5.19|) holds uniformly for all large A"; 
i.e., that the error term 0(|v| 4 ) can be bounded independently of N. For 
fixed N and a finite box, Eq. (j5.19j) with p s /p = 1 always holds for a Bose 
gas with an arbitrary interaction if v is small enough, owing to the discrete- 
ness of the energy spectrum. 1 There are other definitions of the superfluid 
density that may lead to different results PrSv , but this is the one we shall 
use here and shall not dwell on this issue since it is not clear that there is a 
"one-size-fits-all" definition of superfluidity. For instance, in the definition 
we use here the ideal Bose gas is a perfect superfluid in its ground state, 
whereas the definition of Landau in terms of a linear dispersion relation of 
elementary excitations would indicate otherwise. Our main result is that 
with the definition adopted here there is 100% superfluidity in the ground 
state of a 3D Bose gas in the GP limit explained in the previous subsection. 

One of the unresolved issues in the theory of superfluidity is its relation 
to Bose-Einstein condensation (BEC). It has been argued that in general 
neither condition is necessary for the other (c.f., e.g., |H1 lABCUl IKTj ) . but 
in the case considered here, i.e., the GP limit of a 3D gas, we show that 
100% BEC into the constant wave function (in the rest frame) prevails even 
if an external velocity field is imposed. A simple example illustrating the 
fact that BEC is not necessary for superfluidity is the ID hard-core Bose 
gas. This system is well known to have a spectrum like that of an ideal 
Fermi gas |Gij (see also Section 8), and it is easy to see that it is superfluid 
in its ground state in the sense of ()5.19[) . On the other hand, it has no BEC 
jLel IPiStj . The definition of the superfluid velocity as the gradient of the 



The ground state with v = remains an eigenstate of the Hamiltonian with arbitrary 
v (but not necesssarily a ground state) since its total momentum is zero. Its energy is 
^mJVv 2 above the ground state energy for v = 0. Since in a finite box the spectrum of 
the Hamiltonian for arbitrary v is discrete and the energy gap above the ground state is 
bounded away from zero for v small, the ground state for v = is at the same time the 
ground state of the Hamiltonian with v if ^mNv 2 is smaller than the gap. 
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phase of the condensate wave function |HoM[ IBmj is clearly not applicable 
in such cases. 

We consider a Bose gas with the Hamiltonian (|2.1|) in a box K, of side 
length L, assuming periodic boundary conditions in all three coordinate 
directions. Imposing an external velocity field v = (0,0, ±|v|) means that 
the momentum operator p = —iftV is replaced by by p — mv, retaining the 
periodic boundary conditions. The Hamiltonian in the moving frame is thus 

JV 

H 'n = -vY, V L+ E (5.20) 

j=l l<i<j<N 

where Vj )¥ , = Vj + i(0, 0, <p/L) and the dimensionless phase ip is connected 
to the velocity v by 

±|v|Lm 

<P = h ■ (5-21) 

Let Eo(N,a,ip) denote the ground state energy of (|5.20j) with periodic 
boundary conditions. Obviously it is no restriction to consider only the case 
— 7r < ip < it, since Eq is periodic in ip with period 2ir (see Remark 1 below). 
For the ground state of H' N , let 7/v be its one-particle reduced density 
matrix. We are interested in the Gross-Pitaevskii (GP) limit N — > oo with 
Na/L fixed. We also fix the box size L. This means that a should vary 
like 1/N which, as explained in the previous subsection, can be achieved 
by writing v(r) = a~ 2 v\{r /a), where v\ is a fixed potential with scattering 
length 1, while a changes with N. 

Theorem 5.3 (Superfluidity and BEC of homogeneous gas). For 

M — 11 

hm — = ATT^ap + /i— r (5.22) 

in the limit N — ► oo with Na/L and L fixed. Here p = N/L 3 , so ap is fixed 
too. In the same limit, for \tp\ < it, 

lim — 7jy(s,f) = (5.23) 

JV— >CXD JV L 6 

in trace class norm, i.e., lirmv^oo Tr [|7jv/jV— |L _3 / 2 )(L _3 / 2 | | ] = 0. 

Note that, by the definition (fFTT3j) of p s and Eq. (jS2B, Eq. (|^2^|) means 
that p s = p, i.e., there is 100% superfluidity. For cp = 0, Eq. (|5.22|) follows 
from Eq. while (15^3]) for p = is the BEC of Theorem O 2 

Remarks. 1. By a unitary gauge transformation, 

(^)(x lj ...,x A r) = e i ^^/ L ^(xi,...,x iV ) , (5.24) 



2 The convention in Theorem 15.11 where p and Na/L stay fixed, is different from the 
one employed here, where L and Na/L are fixed, but these two conventions are clearly 
equivalent by scaling. 
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the passage from (j'i.lj) to (j5.2()j) is equivalent to replacing periodic boundary 
conditions in a box by the twisted boundary condition 

#0?1 + (0, 0, L), x 2 , . . . , x N ) = e iv ^(x i, x 2 , ...,x N ) (5.25) 

in the direction of the velocity field, while retaining the original Hamiltonian 

J231). 

2. The criterion \p\ < tt means that |v| < irh/(mL). The corresponding 
energy ^m(irH/(mL)) 2 is the gap in the excitation spectrum of the one- 
particle Hamiltonian in the finite-size system. 

3. The reason that we have to restrict ourselves to \p>\ < tt in the second 
part of Theorem 15.31 is that for \tp\ = tt there are two ground states of the 
operator (V + i(f/L) 2 with periodic boundary conditions. All we can say 
in this case is that there is a subsequence of 77V that converges to a density 
matrix of rank < 2, whose range is spanned by these two functions 

Proof of Theorem As in the proof of Theorem 15.11 we combine the lo- 
calization Lemma 15.21 this time with p 7^ 0, and a generalized Poincare 
inequality, this time Lemma 14.21 We thus infer that, for any symmetric ^ 
with V?) = 1 and for TV large enough, 

l(*,i^)(l- const. yVl7)-i 



c 



+-7T / dX dx 
^ Jk."- 1 Jk 



1 



^(xi, X) — L~ 3 [J K dx^(x, X) 



2 



where we used that |J7 C | < ^-NR 3 = const. L 3 Y 2 / 17 . From this we can 
infer two things. First, since the kinetic energy, divided by N, is certainly 
bounded independently of N, as the upper bound shows, we get that 

r E (N,a,p) Lo 2 
lim mf uv ' ,yj > An/ipa + ^ (5.26) 
N-^oo N L 2 

for any \p\ < tt. By continuity this holds also for \p\ = tt, proving (|5.22j) . (To 
be precise, Eq/N — fip 2 L~ 2 is concave in p, and therefore stays concave, and 
in particular continuous, in the limit N — » 00.) Secondly, since the upper 
and the lower bounds to Eq agree in the limit considered, the positive last 
term in (j5.lt) jl has to vanish in the limit. I.e., we get that for the ground 
state wave function v&o of H' N 



lim 

N—fOO 



[ dX I dxx ^ (xi,X) - L- 3 [Ldx^ (x,X)] 2 = . (5.27) 

Using again (j5.18j) . this proves (|5.23|) in a weak sense. As explained in 
|LSe| ILSSYj . this suffices for the convergence N _1 jn —> |L~ 3 / 2 )(L~ 3 / 2 | in 
trace class norm. ■ 
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Theorem 15, 31 can be generalized in various ways to a physically more re- 
alistic setting, for example replacing the periodic box by a cylinder centered 
at the origin. We shall comment on such extensions at the end of Sectional 

6. Gross-Pitaevskii Equation for Trapped Bosons 

In the recent experiments on Bose condensation (see, e.g., |KD| ) . the par- 
ticles are confined at very low temperatures in a 'trap' where the particle 
density is inhomogeneous, contrary to the case of a large 'box', where the 
density is essentially uniform. We model the trap by a slowly varying con- 
fining potential V, with V(x) — > oo as \x\ — * oo. The Hamiltonian becomes 

N 

H = Y J {-^ + V(x i )}+ Yl v(\xi-Xj\). (6.1) 

i=l l<i<j<N 

Shifting the energy scale if necessary we can assume that V is nonnegative. 
The ground state energy, Hu, of — + V(x) is a natural energy unit and 
the corresponding length unit, y hf (mu) = y / 2fj,/(huj) = L osc , is a measure 
of the extension of the trap. 

In the sequel we shall be considering a limit where a/L osc tends to 
zero while N — ► oo. Experimentally a/L osc can be changed in two ways: 
One can either vary Lose o*" The first alternative is usually simpler 
in practice but very recently a direct tuning of the scattering length it- 
self has also been shown to be feasible [CCRCW . Mathematically, both 
alternatives are equivalent, of course. The first corresponds to writing 
V{x) = L~^ c V\{x / L osc ) and keeping V\ and v fixed. The second corresponds 
to writing the interaction potential as f(|x|) = a~ 2 v\(\x\/ a) like in (|5.4|) . 
where v\ has unit scattering length, and keeping V and V\ fixed. This is 
equivalent to the first, since for given V\ and v% the ground state energy of 
(|6.1j) . measured in units of hw, depends only on N and a/L osc . In the dilute 
limit when a is much smaller than the mean particle distance, the energy 
becomes independent of v\. 

We choose L osc as a length unit. The energy unit is hu = 2fiL^ c = 2fi. 
Moreover, we find it convenient to regard V and V\ as fixed. This justifies 
the notion Eq (N, a) for the quantum mechanical ground state energy. 

The idea is now to use the information about the thermodynamic limiting 
energy of the dilute Bose gas in a box to find the ground state energy of 
(|6.1j) in an appropriate limit. This has been done in |LSeYl| |LSeY2j and in 
this section we give an account of this work. As we saw in Sections El and 131 
there is a difference in the p dependence between two and three dimensions, 
so we can expect a related difference now. We discuss 3D first. 

6.1. Three Dimensions. Associated with the quantum mechanical ground 
state energy problem is the Gross-Pitaevskii (GP) energy functional |Grl| 
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EH El 

£ GP [<f>\= [ {v\V(f>\ 2 + V\4>\ 2 + ^ l ia\(f ) \ i )dx (6.2) 
with the subsidiary condition 

\4>\ 2 = N. (6.3) 



As before, a > is the scattering length of v. The corresponding energy is 

E GP (N,a) = inf £ GP [</,] = £ GP [0 GP ], (6.4) 
f\<t>\*=N 

with a unique, positive 4> GP . The existence of the minimizer <j) is proved 
by standard techniques and it can be shown to be continuously differentiable, 
see |LSeYlj . Sect. 2 and Appendix A. The minimizer depends on N and a, 
of course, and when this is important we denote it by (fiWa' 

The variational equation satisfied by the minimizer is the GP equation 

-fiA(j) GP (x) + V(x)<tP p {x) + 8vr/xa(/) GP (x) 3 = p GP > GP (x) , (6.5) 

where fx is the chemical potential, given by 

fi GP = dE GP (N,a)/alN = E GP {N,a)/N + (47r^a/N) J |0 GP (x)| 4 dx. (6.6) 

The GP theory has the following scaling property: 

E GP (N, a) = NE GP (l, Na), (6.7) 

and 

= NV 2 <j> GP Na (x). (6.8) 

Hence we see that the relevant parameter in GP theory is the combination 
Na. 

We now turn to the relation of E GP and (f> GP to the quantum mechanical 
ground state. If v = 0, then the ground state of 1)6. If) is 

V (xi, x N ) = UiLiMxi) 

with (fio the normalized ground state of — + V(x). In this case clearly 
(j) GP = y/~N 4>o, and then E GP = Nftw = Eq. In the other extreme, if 
V(x) = for x inside a large box of volume L 3 and V(x) = oo otherwise, 
then 4> GP w yjN/L 3 and we get £ GP (N,a) = 47r/xaA r2 /L 3 , which is the 
previously considered energy Eq for the homogeneous gas in the low density 
regime. (In this case, the gradient term in £ GP plays no role.) 
In general, we expect that for dilute gases in a suitable limit 

E ^E GP and p QM (x) ^ ^ (x)\ 2 ^ p GP (x) , (6.9) 

where the quantum mechanical particle density in the ground state is defined 
by 

p Q M (x) = iV [ \y (x,X2,...,x N )\ 2 dx 2 ---dx N . (6.10) 
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Dilute means here that 

pa 3 < 1, (6.11) 

where 

P=jjj \p GP m 2 dx (6.12) 

is the mean density. 

The limit in which Q6.9JI can be expected to be true should be chosen so 
that all three terms in £ GP make a contribution. The scaling relations 1)6.7)1 
and 1)6.8)1 indicate that fixing Na as N — > oo is the right thing to do (and this 
is quite relevant since experimentally N can be quite large, 10 6 and more, 
and Na can range from about 1 to 10 4 [DGPSj ). Fixing Na (which we refer 
to as the GP case) also means that we really are dealing with a dilute limit, 
because the mean density p is then of the order N (since p~N,a = ^Pi,Na) 
and hence 

a 3 p~N- 2 . (6.13) 
The precise statement of 1)6.9)1 is: 

Theorem 6.1 (GP limit of the QM ground state energy and den- 
sity). If N — > oo with Na fixed, then 

hm Zi^K = 1, (6-14) 



and 

N-+00 N 

in the weak L 1 -sense. 



J™ ^ T PTaW = \&a(x)\ 2 (6.15) 



Convergence can not only be proved for the ground state energy and 
density, but also for the individual energy components: 

Theorem 6.2 (Asymptotics of the energy components). Let ipo de- 
note the solution to the zero-energy scattering equation for v (under the 
boundary condition limm^^ tpo(x) = I) and s = J \ Vipo\ 2 /(Ana). Then 
< s < 1 and, in the same limit as in Theorem \6.1\ above, 



lim / \V 3l ^Q{xx,X)\ 2 dxidX 

\V^ p Na (x)\ 2 dx + 47rNas I \^ p Na (x)\ A dx, (6.16a) 



lim / V(xx)\V (xi,X)\ 2 dxidX = V(x)\(j)^ a (x)[ 2 dx, (6.16b) 

N 



lim / v(\xi - Xj\)\^ (xi,X)\ 2 dxidX 

(l-s)4vriVa f |</>f^ a (f)| 4 dx. (6.16c) 



N- - - 

i=2 
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Here we introduced again the short hand notation (|5.2j) . Theorem 16. 21 is a 
simple consequence of Theorem 16. II bv variation with respect to the different 
components of the energy, as was also noted in [CS2 . More precisely, Eq. 
.141) can be written as 



lim ±-E (N,a) = E GP (l,Na). (6.17) 

N— >oo iv 

The ground state energy is a concave function of the mass parameter p, so 
it is legitimate to differentiate both sides of (|6.17|) with respect to p. In 
doing so, it has to be noted that Na depends on p through the scattering 
length. Using (|2.13j) one sees that 



d(pa) 



J \Vip \ 2 dx (6.18) 



dp 

by the Feynman-Hellmann principle, since tpo minimizes the left side of 
(23HJ). 

We remark that in the case of a two-dimensional Bose gas, where the 
relevant parameter to be kept fixed in the GP limit is N/\ ln(a 2 p7v)| (c.f. 
Sections IHl and l6~2l ) . the parameter s in Theorem 16.21 can be shown to be 
always equal to 1. I.e., in 2D the interaction energy is purely kinetic in the 
GP limit (see jESU)- 



To describe situations where Na is very large, it is appropriate to consider 
a limit where, as N —* oo, a 3> N~ l , i.e. Na —* oo, but still pa? —* 0. In this 
case, the gradient term in the GP functional becomes negligible compared 
to the other terms and the so-called Thomas-Fermi (TF) functional 

£ TF [p] = [ (Vp + 4Trpap 2 )dx (6.19) 

Jr3 

arises. (Note that this functional has nothing to do with the fermionic 
theory invented by Thomas and Fermi in 1927, except for a certain formal 
analogy.) It is defined for nonnegative functions p on M 3 . Its ground state 
energy E TF and density p TF are defined analogously to the GP case. (The 
TF functional is especially relevant for the two-dimensional Bose gas. There 
a has to decrease exponentially with N in the GP limit, so the TF limit is 
more adequate; see Subsection 16. 21 below) . 

Our second main result of this section is that minimization of (|6,19|) re- 
produces correctly the ground state energy and density of the many-body 
Hamiltonian in the limit when N — » oo, a s p — » 0, but Na — > oo (which we 
refer to as the TF case), provided the external potential is reasonably well 
behaved. We will assume that V is asymptotically equal to some function 
W that is homogeneous of some order s > 0, i.e., W{\x) = X s W(x) for all 
A > 0, and locally Holder continuous (see |LSeY2j for a precise definition). 
This condition can be relaxed, but it seems adequate for most practical 
applications and simplifies things considerably. 
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Theorem 6.3 (TF limit of the QM ground state energy and den- 
sity). Assume that V satisfies the conditions stated above. If g = Na — > oo 
as N — > oo, but still a?p — > 0, i/ien 

j te, ^(jV,a) =1 ' (6 - 20) 



and 



li- ^^P^(5 1/(s+3) f) = (6.21) 



zn £foe weak L l -sense, where pfj is the minimizer of the TF functional under 
the condition f p = 1, a = 1, and with V replaced by W . 

In the following, we will present the essentials of the proofs Theorems 16. II 
and 16,31 We will derive appropriate upper and lower bounds on the ground 
state energy Eq. 

The proof of the lower bound in Theorem 16.11 presented here is a modified 
version of (and partly simpler than) the original proof in [LSeYlj . 

The convergence of the densities follows from the convergence of the en- 
ergies in the usual way by variation with respect to the external potential. 
For simplicity, we set p = 1 in the following. 

Proof of Theorems 16'. il and \6'.'A Part 1: Upper bound to the QM energy. To 
derive an upper bound on Eq we use a generalization of a trial wave function 
of Dyson |Dlj . who used this function to give an upper bound on the ground 
state energy of the homogeneous hard core Bose gas (c.f. Section l2~T|) . It is 
of the form 

N 

V(xx,...,x N ) = l[4> GP (x l )F(x 1 ,...,x N ), (6.22) 

i=l 

where F is constructed in the following way: 

N 

F(xx, x N ) = Y[ f{ti{x\, Xi)), (6.23) 

i=l 

where t{ = min{| |, 1 < j < i — 1} is the distance of Xi to its nearest 

neighbor among the points x±, . . . , Xi—i, and / is a function of r > 0. As in 
(|2.20f) we choose it to be 

f(r) = { /0(r)/ f° (6) £ \l\ (6.24) 



where /o is the solution of the zero energy scattering equation (|2.3|) and b 
is some cut-off parameter of order b ~ p^ 1 ^ 3 . The function (|6.22j) is not 
totally symmetric, but for an upper bound it is nevertheless an acceptable 
test wave function since the bosonic ground state energy is equal to the 
absolute ground state energy. 
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The result of a somewhat lengthy computation (see |LSeYlj for details) 
is the upper bound 

E (N, a) < E GP (N, a) (l + 0{ap 1/3 )^j . (6.25) 



Part 2: Lower bound to the QM energy, GP case. To obtain a lower bound 
for the QM ground state energy the strategy is to divide space into boxes and 
use the estimate on the homogeneous gas, given in Theorem l2.41 in each box 
with Neumann boundary conditions. One then minimizes over all possible 
divisions of the particles among the different boxes. This gives a lower bound 
to the energy because discontinuous wave functions for the quadratic form 
defined by the Hamiltonian are now allowed. We can neglect interactions 
among particles in different boxes because v > 0. Finally, one lets the 
box size tend to zero. However, it is not possible to simply approximate 
V by a constant potential in each box. To see this consider the case of 
noninteracting particles, i.e., v = and hence a = 0. Here Eq = NtkJ, but 
a 'naive' box method gives only mm^V(x) as lower bound, since it clearly 
pays to put all the particles with a constant wave function in the box with 
the lowest value of V. 

For this reason we start by separating out the GP wave function in each 
variable and write a general wave function ^ as 



N 



^(xi 



(6.26) 



Here 



^ p is normalized so that / |^ GP | 2 = N. Eq. (|6~26|) defines F 
for a given because 4> GP is everywhere strictly positive, being the ground 
state of the operator — A + y + 87ra|<^> GP | 2 . We now compute the expectation 
value of H in the state ^f. Using partial integration and the variational 
equation (j6.5|) for 4> , we see that 



with 



(*|fT*) 



E GP (N,a) = 4vra 



„GP|2 



+ Q(F), 



(6.27) 



Q(F) = 

n ml 



E 

i=l 



P GF (x k ) ( \ViF\ 2 + lE^XWi ~ %l) " 8™ P GP (xi) 



/nf=iP GP (^)i^i 



\F\ 



(6.28) 



We recall that p 



GP, 



x 



x 



For computing the ground state energy 



of H we have to minimize the normalized quadratic form Q. Compared to 
the expression for the energy involving ^ itself we have thus obtained the 
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replacements 

N N 

V{x) -8vrap GP (x) and \\dxi -» J^ / o GP (f i )^ • (6-29) 

i=i i=i 

We now use the box method on this problem. More precisely, labeling the 
boxes by an index a, we have 

infQ(F) > inf VinfQ a (F Q ), (6.30) 

F {n a \ — ' Fa 

where Q a is defined by the same formula as Q but with the integrations 
limited to the box a, F a is a wave function with particle number n a , and 
the infimum is taken over all distributions of the particles with ^ n a = N. 

We now fix some M > 0, that will eventually tend to oo, and restrict 
ourselves to boxes inside a cube Am of side length M. Since v > the 
contribution to (|6.30|) of boxes outside this cube is easily estimated from 
below by — 87riVa sup^ AM P GF (%)i w hich, divided by N, is arbitrarily small 
for M large, since Na is fixed and <f> GP /N 1 ^ 2 = GP /„ decreases faster than 
exponentially at infinity ( L SeYlj . Lemma A. 5). 

For the boxes inside the cube Am we want to use Lemma [2.5l and therefore 
we must approximate p GP by constants in each box. Let p ajT aax an d p a ,mm, 
respectively, denote the maximal and minimal values of p GP in box a. Define 

^a(xi,. . . ,f„J = F a (xi, . . .,f„J Y[ 0° P (^fc)) ( 6 - 31 ) 



and 



k=l 



^\x 1 ,...,x n J=F a (x 1 ,...,x n j]Jcp GP (x k ). (6.32) 



k=l 

k^i 



(6.33) 



We have, for all 1 < % < n a , 

J II P GP (Zk) (|V 4 F Q | 2 + v(\xi - x-MFaf) 

k=l j^i 

> Pa*** J (iv^»i 2 + \ j2 v(\3t - m^f) ■ 

We now use Lemma 12.51 to get, for all < e < 1, 

$M>P*MnJ (e\Vi*®\ 2 + a(l-e)U(U)\*^ 2 ) (6.34) 

where U is the distance to the nearest neighbor of Xi, c.f., I|2.29j) . and U the 
potential (l2~3oT) . 

Since ^ a = (p GP (x^i&a we can estimate 

|V^ Q | 2 < 2 Pa ^ x \vM: ] \ 2 + 2\^\ 2 NC M (6.35) 
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with 

Cm = ^ sup |V0 GP (x)| 2 = sup |V^P a (f)| 2 . (6.36) 

JV xeAjif x&A M 
Since Na is fixed, Cm is independent of N. Inserting (|6.35f) into (|6.34|) . 
summing over i and using p GP (x{) < p a , m ax in the last term of (|6,28|) (in 
the box a), we get 

Qa{Fa) > Pa ' mm E^(n a , L) - 8irap aima , x n a - sCun a , (6.37) 

Pa, max 

where L is the side length of the box and E^(n a ,L) is the ground state 
energy of 

^(-±eAi + (l-e)aU(U)) (6.38) 

i=l 

in the box (c.f. (|2.35|) ). We want to minimize (|6.37|) with respect to n a and 
drop the subsidiary condition Y^ a n a = N in (|6.30|) . This can only lower the 
minimum. For the time being we also ignore the last term in (|6.37|) . (The 
total contribution of this term for all boxes is bounded by eCmN and will 
be shown to be negligible compared to the other terms.) 

Since the lower bound for the energy of Theorem 12.41 was obtained pre- 
cisely from a lower bound to the operator 1)6.38(1 . we can use the statement 
and proof of Theorem 12.41 From this we see that 

E^(n a ,L) > (1 - e)^(l - CYW) (6.39) 

with Y a = a 3 n Q /L 3 , provided Y a is small enough, e > Y<t /17 and 
( const. )Y a . The condition on e is certainly fulfilled if we choose e = 
y 1 / 17 with Y = a?N/L 3 . We now want to show that the n a minimizing the 
right side of Q6.37JI is large enough for (|6.39|) to apply. 

If the minimum of the right side of (|6.37|) (without the last term) is taken 
for some n a , we have 



Pa,m'm / „u 



Pa, max 

On the other hand, we claim that 
Lemma 6.4. For any n 



E v £ {n a + 1, L) - E?(n a ,L)) > 8vrap a , max . (6.40) 



E^(n + 1,L) - E^{n,L) <8na^. (6.41) 

Proof. Denote the operator (|6.38f) by H n , with n a = n, and let ^ n be its 
ground state. Let be the distance to the nearest neighbor of Xj among 
the n + 1 points x%, . . . , x n+ \ (without Xi) and U the corresponding distance 
excluding x n+ \. Obviously, for 1 < i < n, 

U{ft<U(t i ) + U{\a! i -a! n +i\) (6.42) 
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and 

n 

U{t' n+1 )<Y,U{\xi-x n+1 \). (6.43) 



i=l 



Therefore 



H n+ i < H n - ±eA n+1 + 2aJ2 u (Wi ~ *n+i\)- (6.44) 

i=i 

Using ^ n /L 3//2 as trial function for H n +i we arrive at (|6.41|) . ■ 
Eq. (|6.41|) together with (j6.40|) shows that n a is at least ~ Pa,max 

L 3 . We 

shall choose L ~ A^ 1 / 10 , so the conditions needed for (|6.39l) are fulfilled for 
N large enough, since p a)max ~ N and hence n a ~ N 7 / 10 and Y a ~ A^~ 2 . 
In order to obtain a lower bound on Q a we therefore have to minimize 

4to (*H!*L2| (l _ CYW) - 2n aPa , max ) . (6.45) 

We can drop the requirement that n a has to be an integer. The minimum 
of (|6.45|) is obtained for 

o 2 T 3 

na= Pm (i-cY^y (6 - 46) 

By Eq. (|6.27f) this gives the following lower bound, including now the last 
term in (j6.37j) as well as the contributions from the boxes outside Am, 

E (N,a) - E GP (N,a) > 

4™ f \p GP \ 2 - 4™ Y, PlM^ 1 ) 

J aCA M \^Q,min V 1 ^ J / 

- Y 1/17 NC M ~ ^naN sup p GP {x). 

x<£A M 

Now /? GP is differentiable and strictly positive. Since all the boxes are in the 
fixed cube Am there are constants C < oo, C" > 0, such that 

Ptt,max - Pa,mm < NC'L, p Qjm i n > NC". (6.48) 

Since L ~ A^ 1 / 10 and Y ~ A^ 17 / 10 we therefore have, for large N, 
Also, 

4™ £ ^ min L 3 < 4vra / |p GP | 2 < E GP (N,a). (6.50) 

oCA M 

Hence, noting that E GP (N, a) = NE GP (l, Na) ~ N since Na is fixed, 
Ea(N,a) > 1 _ (congt + ^^-1/10 _ ( const .) sup ^J 2 , (6.51) 



E GP (N, a) ~ ' — ^ ; S Ja m 
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where the constants depend on Na. We can now take N — > oo and then 
M -> oo. 

Pari 3: Lower bound to the QM energy, TF case. In the above proof of 
the lower bound in the GP case we did not attempt to keep track of the 
dependence of the constants on Na. In the TF case Na — > oo, so one would 
need to take a closer look at this dependence if one wanted to carry the 
proof directly over to this case. But we don't have to do so, because there is 
a simpler direct proof. Using the explicit form of the TF minimizer, namely 

p T N F a (x) = ^[p TF -V(x)] + , (6.52) 

where = max{t, 0} and p TF is chosen so that the normalization condi- 
tion J* pj^ a = N holds, we can use 

V(x) > p TF - 87rap TF (x) (6.53) 

to get a replacement as in (16.29(1 . but without changing the measure. More- 
over, p TF has compact support, so, applying again the box method described 
above, the boxes far out do not contribute to the energy. However, /i TF 
(which depends only on the combination Na) tends to infinity as Na — ► oo. 
We need to control the asymptotic behavior of /i TF , and this leads to the 
restrictions on V described in the paragraph preceding Theorem 16.31 For 
simplicity, we shall here only consider the case when V itself is homogeneous, 
i.e., V(Xx) = X s V(x) for all A > with some s > 0. 
In the same way as in 1)6.6(1 we have, with g = Na, 

p TF (g) = dE TF (N,a)/dN = E TF (l,g) + ^g J \pJ F g (x)\ 2 dx. (6.54) 

The TF energy, chemical potential and minimizer satisfy the scaling relations 

£ TF (l, ff )=^+ 3 )£ TF (l,l), (6.55) 

^ TF (5)=5 S/{S+3 V TF (1), (6-56) 

and 

g ^pJ F (gV^x) = P J F (x). (6.57) 
We also introduce the scaled interaction potential, v, by 

v(x) = g 2 ^ s+ ^v(g 1/{s+3) x) (6.58) 

with scattering length 

a = g - 1 /^a. (6.59) 
Using ()6.53() . ()6.54() and the scaling relations we obtain 

E (N, a) > E TF (N, a) + ^Ng s '^ f \pJ F \ 2 + g~ 2 '^Q (6.60) 



42 E.H. LIEB, R. SEIRINGER, J. P. SOLOVEJ, AND J. YNGVASON 

with 

vyjbi — Xj)\^\ — on i \ api i \xi ) \ |2 1 



Q= mf_ [ (|Vi*| 2 +iX;^-^)l*l a - ^^2(^)1*1 



(6-61) 

We can now proceed exactly as in Part 2 to arrive at the analogy of Eq. 
(|6.47[) . which in the present case becomes 

E (N,a)-E TF (N,a) > 

4xNg*/W [ \ P J F \ 2 - teNa^fi^O. - CY^y\ (6 - 62) 

Here p a ,ma,x is the maximum of pJ F in the box a, and Y = a 3 N/L 3 . This 
holds as long as L does not decrease too fast with N. In particular, if L is 
simply fixed, this holds for all large enough N. Note that 

p = Np 1>g ~Ng- 3 tt s+ Vp 1>1 , (6.63) 

so that a 3 N ~ a?p goes to zero as N — > oo by assumption. Hence, if we first 
let N — > oo (which implies Y — > 0) and then take L to zero, we of arrive at 
the desired result 

liminf E ^ N > a) >1 (6.64) 

in the limit iV — > oo, a 3 /) — > 0. Here we used the fact that (because V, and 
hence /o TF , is continuous by assumption) the Riemann sum J2 a p 2 a mSbX L 3 
converges to f \pf F \ 2 as L — > 0. Together with the upper bound Q6.25JI 
and the fact that E GP (N , a) / E TF (N , a) = E GP {1, Na)/ E TF (1, Na) -» 1 as 
iVa — » oo, which holds under our regularity assumption on V (c.f. Lemma 
2.3 in |LSeY2| ). this proves (l6~Ti]) and (l6~2T7|) . 



Part ^: Convergence of the densities. The convergence of the energies im- 
plies the convergence of the densities in the usual way by variation of the 
external potential. We show here the TF case, the GP case goes analogously. 
Set again g = Na. Making the replacement 

V(x) — ► V(x) + 5g sl{s+3) Z{g- l ^ s+ ^x) (6.65) 

for some positive Z £ and redoing the upper and lower bounds we see 
that 1)6. 20j) holds with W replaced by W + SZ. Differentiating with respect 
to 5 at 5 = yields 

lim fig^ = pJJW (6-66) 

in the sense of distributions. Since the functions all have L 1 -norm 1, we can 
conclude that there is even weak L 1 -convergence. ■ 
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6.2. Two Dimensions. In contrast to the three-dimensional case the en- 
ergy per particle for a dilute gas in two dimensions is nonlinear in p. In view 
of Schick's formula (|3.1|) for the energy of the homogeneous gas it would ap- 
pear natural to take the interaction into account in two dimensional GP 
theory by a term 

4tt [ llnd^lV)!- 1 !^! 4 ^, (6.67) 
and such a term has, indeed, been suggested in |Shj and KNSQ . However, 



since the nonlinearity appears only in a logarithm, this term is unnecessarily 
complicated as far as leading order computations are concerned. For dilute 
gases it turns out to be sufficient, to leading order, to use an interaction 
term of the same form as in the three-dimensional case, i.e, define the GP 
functional as (for simplicity we put /i = 1 in this section) 

£ GP M= / (|V(/>| 2 + y|(/>| 2 + 47ra|0| 4 )(if, (6.68) 
Jr. 2 

where instead of a the coupling constant is now 

a= llnQ^va 2 )! -1 ( 6 - 69 ) 
with pN the mean density for the GP functional at coupling constant 1 and 
particle number N . This is defined analogously to (|6.12|) as 

PN = l I \&dx (6.70) 



N _ 

where ^ is the minimizer of (|6.68|) with a = 1 and subsidiary condition 
f\(f)\ 2 = N. Note that a in (|6.69j) depends on N through the mean density. 

Let us denote the GP energy for a given N and coupling constant a by 
E GP (N, a) and the corresponding minimizer by <^y^- As in three dimensions 
the scaling relations 

E GP (N,a) = NE GP (l,Na) (6.71) 

and 

N-^ G N P a = ffi al (6.72) 
hold, and the relevant parameter is 

g = Na. (6.73) 

In three dimensions, where a = a, it is natural to consider the limit iV — ► 
oo with g = Na= const. The analogue of Theorem l6.1l in two dimensions is 

Theorem 6.5 (Two-dimensional GP limit theorem). If, for N — > oo ; 

a 2 p~N with g = N/\ ln(a 2 pAr)| fixed, then 

N™ooE GP (N,l/\ln(a 2 p N )\) { ' 

and 

lim lp^(f) = |^(x)| 2 (6.75) 
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in the weak L 1 -sense. 

This result, however, is of rather limited use in practice. The reason is 
that in two dimensions the scattering length has to decrease exponentially 
with N if g is fixed. The parameter g is typically very large in two dimensions 
so it is more appropriate to consider the limit N — > oo and g — > oo (but still 
p N a 2 -»■ 0). 

For potentials V that are homogeneous functions of x, i.e., 

V(Xx) = X s V(x) (6.76) 

for some s > 0, this limit can be described by the a 'Thomas-Fermi' energy 
functional like ()6.19)) with coupling constant unity: 

£ TF [p] = [ (V(x)p(x) + 4irp(x) 2 )dx. (6.77) 

This is just the GP functional without the gradient term and a = 1. Here 
p is a nonnegative function on R 2 and the normalization condition is 

= 1. (6.78) 

The minimizer of (|6.77|) can be given explicitly. It is 

P J f 1 (x) = (8tt)- 1 [p tf -V(x)] + (6.79) 

where the chemical potential p TF is determined by the normalization con- 
dition 1)6.78)) and [t] + = t for t > and zero otherwise. We denote the 
corresponding energy by E TF (l,l). By scaling one obtains 

lim E GP (l,g)/g s / (s+ V = E TF (1, 1), (6.80) 



lim g 2 ^ P r g {g ll{s+2) *) = p3(*). (6-81) 

g— >oo ,a ' 

with the latter limit in the strong L 2 sense. 

Our main result about two-dimensional Bose gases in external potentials 
satisfying 1)6.76 J) is that analogous limits also hold for the many-particle 
quantum mechanical ground state at low densities: 

Theorem 6.6 (Two-dimensional TF limit theorem). In two dimen- 
sions, if a 2 p~N — ► 0, but g = N/\ ln(pAra 2 )| —> oo as N — > oo then 

Hm Eo{N^ = eTF 

and, in the weak L 1 sense, 

lim g ^^ P Ta(9 1/{S+2) x) = P%1&. (6-83) 

l V — yoo iV 
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Remarks: 1 . As in Theorem 16.31 it is sufficient that V is asymptotically 
equal to some homogeneous potential, W. In this case, E TF (l,l) and pj\ 
in Theorem 16.61 should be replaced by the corresponding quantities for W. 
2. From Eq. (|6~ST1) it follows that 

p N ~ iV s /( s + 2 ) (6.84) 

for large N. Hence the low density criterion a 2 p<^ 1, means that a/L osc <C 

N -s/2{s+2)_ 

We shall now comment briefly on the proofs of Theorems 16.51 and 16.61 
mainly pointing out the differences from the 3D case considered previously. 

The upper bounds for the energy are obtained exactly in a same way as in 
three dimensions. For the lower bound in Theorem 16.51 the point to notice 
is that the expression (|6.45|) . that has to be minimized over n a , is in 2D 
replaced by 

I Pa, min f^cx 1 ( y ^ A a Pa, max 



p a , max L 2 |ln(a 2 n a /L 2 )| V I Ha 2 N/L^ J \ln(a 2 p N ) 

(6.85) 

since Eq. (j6.39j) has to be replaced by the analogous inequality for 2D (c.f. 
(|3.31[) ). To minimize (|6.85|) we use the following lemma: 

Lemma 6.7. For < x, b < 1 and k > 1 we have 



x 



2 



b , b 2 ( 1 



2— -xk > -— I 1 + — — - J k 2 . (6. 



\lnx\ \lnb\ ~ \lnb\ \ (2|ln6|) 

Proof. Replacing x by xk and using the monotonicity of In we see that it 
suffices to consider k = 1. Since In a; > —^x~ d for all d > we have 

phibl _ 2 x \_^b[ edx 2 +d _ 2x > c{d){b d edllnbl) -i/(i + d) (6 8?) 
cr \ mx\ b b z b 

with 

c(d) = 2( 2 + d )/( 1+d ) ( - - ^ > -1 - -d 2 

{ > V(2 + d)( 2 + rf )/a+ rf ) (2 + d)V(i+<v - 



Choosing d = 1/| lnb| gives the desired result. 

Applying this lemma with x = a 2 n a /L 2 , b = a 2 /9 Q)max and 



(6i 



, _ Pa,max A C \ | ln(a 2 /9 a;max ) | 

"pa.nrinV |ln(a 2 iV/L 2 )|V5; |ln(a 2 ^)| l ' j 

we get the bound 

W > -47r .^; m 2 axL ' f 1 + — - -J J) k. (6.90) 

|ln(a 2 piv)| V 4|ln(a 2 p Qjmax )| 2 y 

In the limit considered, k and the factor in parenthesis both tend to 1 and 
the Riemann sum over the boxes a converges to the integral as L — > 0. 
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The TF case, Thm. lfTHl is treated in the same way as in three dimensions, 
with modifications analogous to those just discussed when passing from 3D 
to 2D in GP theory. 

7. Bose-Einstein Condensation and Superfluidity for Dilute 

Trapped Gases 

It was shown in the previous section that, for each fixed Na, the mini- 
mization of the GP functional correctly reproduces the large N asymptotics 
of the ground state energy and density of H - but no assertion about BEC 
in this limit was made. We will now extend this result by showing that 
in the Gross-Pitaevskii limit there is indeed 100% Bose condensation in the 
ground state. This is a generalization of the homogeneous case considered in 
Theorem 15 . 1 1 and although it is not the same as BEC in the thermodynamic 
limit it is quite relevant for the actual experiments with Bose gases in traps. 
In the following, we concentrate on the 3D case, but analogous considera- 
tions apply also to the 2D case. We also discuss briefly some extensions of 
Theorem 15.31 pertaining to superfluidity in trapped gases. 

As in the last section we choose to keep the length scale L osc of the con- 
fining potential fixed and thus write Na instead of Na/L osc . Consequently 
the powers of TV" appearing in the proofs are different from those in the proof 
Theorem 15.11 where we kept Na/L and iV/L 3 fixed. 

For use later, we define the projector 

pGP^GP^GPj (71) 

Here (and everywhere else in this section) we denote 4> GP = G y\r a for simplic- 
ity, where ^"p^a * s ^ ne minimizer of the GP functional Q6.2JI with parameter 
Na and normalization condition J \ cp\ 2 = 1 (compare with (|6.8j0 . Moreover, 
we set /i = 1. 

In the following, denotes the (nonnegative and normalized) ground 
state of the Hamiltonian (|6,1|) . BEC refers to the reduced one-particle den- 
sity matrix j(x, x 1 ) of *S>q, defined in ()5.1|) . The precise definition of BEC 
is is that for some c > this integral operator has for all large N an an 
eigenf unction with eigenvalue > cN. 

Complete (or 100%) BEC is defined to be the property that ^j(x, x 1 ) not 
only has an eigenvalue of order one, as in the general case of an incomplete 
BEC, but in the limit it has only one nonzero eigenvalue (namely 1). Thus, 
jjl(x, x') becomes a simple product <p(x)* f{x f ) as N — ► oo, in which case 
(p is called the condensate wave function. In the GP limit, i.e., N — ► oo 
with Na fixed, we can show that this is the case, and the condensate wave 
function is, in fact, the GP minimizer c/> GP . 

Theorem 7.1 (Bose-Einstein condensation in a trap). For each fixed 
Na 

Jim ^-7(x,f)=0 GP (x)0 GP (x / ) . 
N— >oo N 
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in trace norm, i.e., Tr — P GP \ — ► 0. 

We remark that Theorem 17.11 implies that there is also 100% condensation 
for all ra-particle reduced density matrices 

(n) / - -> . \ 

T v^-l) • • • j -£n> 2-1; • • • > 

= n! (r^) / *o(»lj--- J »Jv)*o(^l)--->^,«n+lj---^)^n+l- ,, <^ 

(7.2) 

of \&o, i- e -; they converge, after division by the normalization factor, to the 
one-dimensional projector onto the n-fold tensor product of <p GP . In other 
words, for n fixed particles the probability of finding them all in the same 
state (f) tends to 1 in the limit considered. To see this, let a*, a denote the 
boson creation and annihilation operators for the state ^> GP , and observe 
that 

1> lim N~ n Nf \(a*) n a n \^ Q ) = lim N- n {^ \(a*a) n \^ ) , (7.3) 

since the terms coming from the commutators [a, a*] = 1 are of lower order 
as iV — > oo and vanish in the limit. From convexity it follows that 

iV- n (^oKa*a) n |*o> > A^ n (^oKa|^o} n (7-4) 

which converges to 1 as N — > oo, proving our claim. 

Another corollary, important for the interpretation of experiments, con- 
cerns the momentum distribution of the ground state. 

Corollary 7.2 (Convergence of momentum distribution). Let 

p(k) = J J j(x, a?) exp[ik ■ (x — x')]dxdx' 
denote the one-particle momentum density o/^o- Then, for fixed Na, 

Jim ^P(k) = $ GP (k)\ 2 

strongly in L 1 (M 3 ). Here, GP denotes the Fourier transform of (j) GP . 

Proof. If T denotes the (unitary) operator 'Fourier transform' and if h is an 
arbitrary L°°-function, then 



= |Tr[^- 1 (7/iV-P GP )^/t]| 
< WhW^Tr^/N-P^], 



from which we conclude that 

\\p/N - |^ GP | 2 ||i < Tr| 7 /iV-P GP | 
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As already stated, Theorem 17.11 is a generalization of Theorem 15. 1| the 
latter corresponding to the case that V is a box potential. It should be noted, 
however, that we use different scaling conventions in these two theorems: In 
Theorem 15.11 the box size grows as iV 1 / 3 to keep the density fixed, while 
in Theorem 17.11 we choose to keep the confining external potential fixed. 
Both conventions are equivalent, of course, c.f. the remarks in the second 
paragraph of Sectional but when comparing the exponents of ./V that appear 
in the proofs of the two theorems the different conventions should be born 
in mind. 

As in Theorem 15.11 there are two essential components of our proof of 
Theorem 17.11 The first is a proof that the part of the kinetic energy that 
is associated with the interaction v (namely, the second term in l|fi.loa|) ) 
is mostly located in small balls surrounding each particle. More precisely, 
these balls can be taken to have radius roughly A r ~ 5 / 9 , which is much smaller 
than the mean-particle spacing A r_1//3 . (The exponents differ from those 
of Lemma 15.21 because of different scaling conventions.) This allows us to 
conclude that the function of x defined for each fixed value of X by 

fx^) = ^-fo(x,X)>0 (7.5) 

has the property that Vgfg(x) is almost zero outside the small balls centered 
at points of X. 

The complement of the small balls has a large volume but it can be a weird 
set; it need not even be connected. Therefore, the smallness of Vgfjt(x) in 
this set does not guarantee that fx{x) is nearly constant (in x), or even that 
it is continuous. We need fx(x) to be nearly constant in order to conclude 
BEC. What saves the day is the knowledge that the total kinetic energy of 
fx(x) (including the balls) is not huge. The result that allows us to combine 
these two pieces of information in order to deduce the almost constancy of 
fjf(x) is the generalized Poincare inequality in Lemma 14.31 The important 
point in this lemma is that there is no restriction on concerning regularity 
or connectivity. 

Using the results of Theorem l6.21 partial integration and the GP equation 
(i.e., the variational equation for </> GP , see Eq. (|6.5[) ) we see that 

lim f \4> GV {x)\ 2 \Vgf^{x)\ 2 dxdX = iirNas [ \<f) GP (x)\ 4 dx . (7.6) 

N^oo J A J 

The following Lemma shows that to leading order all the energy in (|7.fi|) is 
concentrated in small balls. 

Lemma 7.3 (Localization of the energy in a trap). For fixed X let 



= <, x e 



min \x - x k \ > N" 1 ^ 6 } (7.7) 

k>2 1 
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lim I dX I ds|0 GP (£)| 2 |Va/_£(£)| 2 = . 



Remark. In the proof of Theorem l5.1l we chose (5 to be 4/51, but the following 
proof shows that one can extend the range of 5 beyond this value. 



Proof. We shall show that 

dX [ dx |0 GP (x)| 2 |V^(x) 



x 



+ I dX I dx\0 GP (x)\ 2 \U(x)\ 2 



i5>(|x-x fc |)-8^AM0 GP (x)| 2 

k>2 



>-4irNa I \<p GP (x)\ 4 dx-o(l) (7. 



as N — > oo. We claim that this implies the assertion of the Lemma. To see 
this, note that the left side of (|7.8|) can be written as 

l-Eo-fj^-JdxJ dx|0 GP (x)| 2 |V s -/x(^)| 2 , (7-9) 



X 



where we used partial integration and the GP equation (|6.5|) . and also the 
symmetry of \EV The convergence of the energy in Theorem 16.11 and the 
relation (|6,6|) now imply the desired result. 

The proof of (|7.8|) is actually just a detailed examination of the lower 
bounds to the energy derived in |LSeYlj and |LYlj and described in Sec- 
tions and H3 We use the same methods as there, just describing the differ- 
ences from the case considered here. 

Writing 

ff(x) = H<p GP (x k )F(x,X) (7.10) 

k>2 

and using that F is symmetric in the particle coordinates, we see that (|7.8(1 
is equivalent to 



jfQs{F) > -AnNa J |0 GP | 4 - o(l), (7.11) 
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where Q$ is the quadratic form 

AT „ N 




+ E f<\ 



N 



xiDiFi 2 ni^ Gp (^)i 2< ^* 



l<i<j<N J 



k=l 



N 



8nNaY^ / \<t> 



GP 



(x 



^)| 2 |F| 2 ni0 GP (^)| 2 ^fc- (7-12) 



k=l 



Here denotes the set 



= G R 3N \ mm\xi - x k \ < A^ 1 / 3 -* 5 }. 



k+i 



While Q7.11J1 is not true for all conceivable F's satisfying the normalization 
condition 



° k=l 

it is true for an F, such as ours, that has bounded kinetic energy (J7.6|) . 
Looking at Sectional we see that Eqs. (jfHTj) - (joT2Hi (jpTT7)) - (|H3T|) are sim- 
ilar to ()7.11|) . (|7.12|) and almost establish (|7.11|) . but there are differences 
which we now explain. 

In our case, the kinetic energy of particle i is restricted to the subset of 
M 3N in which min^^j \xi — x k \ < N^ 1 ^^ 5 . However, looking at the proof 
of the lower bound to the ground state energy of a homogeneous Bose gas 
discussed in Section 2, which enters the proof of Theorem I6.1[ we see that 
if we choose 5 < 4/51 only this part of the kinetic energy is needed for the 
lower bound, except for some part with a relative magnitude of the order 
e = 0(N~ 2a ) with a = 1/17. (Here we use the a priori knowledge that 
the kinetic energy is bounded by (|7.6|1 .) We can even do better and choose 
some 4/51 < 5 < 2/9, if a is chosen small enough. (To be precise, we choose 
(3 = 1/3 + a and 7 = 1/3 — 4a in the notation of (|2.56|) , and a small enough) . 
The choice of a only affects the magnitude of the error term, however, which 
is still o(l) as N — > 00. ■ 



Proof of Theorem\7l\ For some R > let K = {x G M 3 , \x \ < R}, and 



We shall use Lemma IP1 with d = 3, h(x) = \(p GP (x)\ 2 / J K \c/) GP \ 2 , O = 
Q^nfC and f(x) = fg(x) - (/^}/c (see (J72J) and dZSJ). Since <p GP is 
bounded on K, above and below by some positive constants, this Lemma 
also holds (with a different constant C) with dx replaced by \(f) GP (x)\ 2 dx in 




v 



define 
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Therefore, 



dX 

< C 



GP (^)| 2 |V^(x)| 2 df 



K. 

dX 



R 2 JK 



\cf>^(x)\ 2 \VMx)\ 2 dx 



(7.13) 



where we used that \Q% D JC\ < (47r/3)iV . The first integral on the right 
side of 1)7.13(1 tends to zero as N — > oo by Lemma 17.31 and the second is 
bounded by (|7,6|) . We conclude, since 



l0 GP (^)| 2 /x(^)^< 



because of the positivity of f-g, that 



|<^ GP (x)| 2 /x(^)^ 



1 



liminf-<<^| 7 |</> 



GP\ 



> 



/> GP (£)| 2 dx lim 

iV— »oo 



|^ GP (^)| 2 dx 



dX I dx\V (x,X)\' 

K. 



where the last equality follows from (|6.15j) . Since the radius of K, was arbi- 
trary, we conclude that 

1 



hm -(^P^^GPj 
N^oo TV 



1, 



implying convergence of 7/iV to P GP in Hilbert-Schmidt norm. Since the 
traces are equal, convergence even holds in trace norm (cf. jHU, Thm. 2.20), 
and Theorem 17. II is proven. ■ 

We remark that the method presented here also works in the case of a two- 
dimensional Bose gas. The relevant parameter to be kept fixed in the GP 
limit is N/\ ln(a 2 Pn)\, ah other considerations carry over without essential 
change, using the results in |LSeY2[ lLY2| . cf. Sections El and E21 It should 
be noted that the existence of BEC in the ground state in 2D is not in 
conflict with its absence at positive temperatures |Ho[ IMj . In the hard core 
lattice gas at half filling precisely this phenomenon occurs |KLSj . 

of 

Finally, we remark on generalizations of Theorem 15.31 on superfluidity 
from a torus to some physically more realistic settings. As an example, let 
C be a finite cylinder based on an annulus centered at the origin. Given a 
bounded, real function a(r, z) let A be the vector field (in polar coordinates) 
A(r, 6, z) = ipa(r, z)eg, where eo is the unit vector in the 9 direction. We also 
allow for a bounded external potential V(r, z) that does not depend on 6. 
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Using the methods of Appendix A in LSeYl, it is not difficult to see 
that there exists a p>o > 0, depending only on C and a(r, z), such that for all 
| p | < ipo there is a unique minimizer </> GP of the Gross-Pitaevskii functional 



£ GP [0] = J (|(V + i^(x)) ( / ) (x)| 2 + y(f)|^(f)| 2 + 47r^a|^(f)| 4 )df (7.14) 

under the normalization condition J|</>| 2 = 1. This minimizer does not de- 
pend on 6, and can be chosen to be positive, for the following reason: The 
relevant term in the kinetic energy is T = —r~ 2 [d/d9 + ipr a(r, z)] 2 . If 
\pra(r,z)\ < 1/2, it is easy to see that T > p 2 a(r,z) 2 , in which case, with- 
out raising the energy, we can replace 4> by the square root of the ^-average 
of \4>\ 2 - This can only lower the kinetic energy |LLoj and, by convexity of 
x — ► x 2 , this also lowers the 4 term. 

We denote the ground state energy of £ by E GP , depending on Na and 
p. The following Theorem 17.41 concerns the ground state energy Eq of 

N 

J2[-{V j + iA(x j )) 2 + V(x j )]+ Yl ( 7 - 15 ) 

j=l l<i<j<N 

with Neumann boundary conditions on C, and the one-particle reduced den- 
sity matrix 77V of the ground state, respectively. Different boundary condi- 
tions can be treated in the same manner, if they are also used in (|7.14[) . 

Remark. As a special case, consider a uniformly rotating system. In this 
case A(x) = (preg, where 2ip is the angular velocity. is the Hamiltonian 
in the rotating frame, but with external potential V(x) + A(x) 2 (see e.g. 
[Bm1 p. 131]). 

Theorem 7.4 (Superfluidity in a cylinder). For \ip\ < po 

MN_a lA = EG , {NaM 

N^oo iv 

in the limit N — > 00 with Na fixed. In the same limit, 

ton -UjvO?,*) = GP (f)0 GP (x') (7.17) 

N— >oo iv 

in trace class norm, i.e., limTv^ooTr [|7at/A^— |0 gp )(0 gp | 1=0. 

In the case of a uniformly rotating system, where 2ip is the angular ve- 
locity, the condition \p\ < ip$ in particular means that the angular velocity 
is smaller than the critical velocity for creating vortices |Se31 IFSj . 



Remark. In the special case of the curl-free vector potential A(r, 6) = 
pr~ x e.Q, i.e., a{r,z) = r -1 , one can say more about the role of <po. In 
this case, there is a unique GP minimizer for all tp Z + |, whereas there 
are two minimizers for p G Z + ^ . Part two of Theorem 17.41 holds in this 
special case for all p Z + | , and (|7.16|) is true even for all p. 
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8. One-Dimensional Behavior of Dilute Bose Gases in Traps 

Recently it has become possible to do experiments in highly elongated 
traps on ultra-cold Bose gases that are effectively one-dimensional BBD ( 
IGo| IU1 [S3 IMSKEj . These experiments show peculiar features predicted 
by a model of a one-dimensional Bose gas with repulsive <5-function pair 
interaction, analyzed long ago by Lieb and Liniger |LLj . These include quasi- 
fermionic behavior |Gi| . the absence of Bose-Einstein condensation (BEC) 
in a dilute limit |Le| IPiSt| IGWTj . and an excitation spectrum different 
from that predicted by Bogoliubov's theory |LLl IJKl IKPj . The theoretical 
work on the dimensional cross-over for the ground state in elongated traps 
has so far been based either on variational calculations, starting from a 3D 
delta-potential (OH IDGW1 IGWj . or on numerical Quantum Monte Carlo 
studies |B11 lAGj with more realistic, genuine 3D potentials, but particle 
numbers limited to the order of 100. This work is important and has led 
to valuable insights, in particular about different parameter regions PSW ( 
IDLOj . but a more thorough theoretical understanding is clearly desirable 
since this is not a simple problem. In fact, it is evident that for a potential 
with a hard core the true 3D wave functions do not approximately factorize 
in the longitudinal and transverse variables (otherwise the energy would 
be infinite) and the effective ID potential can not be obtained by simply 
integrating out the transverse variables of the 3D potential (that would 
immediately create an impenetrable barrier in ID). It is important to be 
able to demonstrate rigorously, and therefore unambiguously, that the ID 
behavior really follows from the fundamental Schrddinger equation. It is 
also important to delineate, as we do here, precisely what can be seen in the 
different parameter regions. The full proofs of our assertions are long and 
are given in |LSeY6j . Here we state our main results and outline the basic 
ideas for the proofs. 

We start by describing the setting more precisely. It is convenient to write 
the Hamiltonian in the following way (in units where h = 2m = 1): 

N 

H N)Ljria = ^2(-A j + V r ± (xf) + V L (z j )^+ Yl "ad^-^D 

j=l l<i<j<N 
with x = (x,y, z) = (x , z) and with 

V L {z) = -^V(z/L) , v a (\x\) = ±v(\x\/a) . (8.2) 

Here, r, L, a are variable scaling parameters while V ± , V and v are fixed. 

We shall be concerned with the ground state of this Hamiltonian for large 
particle number N, which is appropriate for the consideration of actual ex- 
periments. The other parameters of the problem are the scattering length, 
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a, of the two-body interaction potential, v, and two lengths, r and L, de- 
scribing the transverse and the longitudinal extension of the trap potential, 
respectively. 

The interaction potential v is supposed to be nonnegative, of finite range 
and have scattering length 1; the scaled potential v a then has scattering 
length a. The external trap potentials V and V ± confine the motion in the 
longitudinal (z) and the transversal (or 1 ) directions, respectively, and are 
assumed to be continuous and tend to oo as \z\ and tend to oo. To 
simplify the discussion we find it also convenient to assume that V is homo- 
geneous of some order s > 0, namely V(z) = \z\ s , but weaker assumptions, 
e.g. asymptotic homogeneity (cf. Section EJ), would in fact suffice. The case 
of a simple box with hard walls is realized by taking s = oo, while the usual 
harmonic approximation is s = 2. It is understood that the lengths associ- 
ated with the ground states of — d 2 /dz 2 + V{z) and — A 1 - + V ± (x ± ) are both 
of the order 1 so that L and r measure, respectively, the longitudinal and 
the transverse extensions of the trap. We denote the ground state energy of 
Q8.1J1 by E® M (N, L, r, a) and the ground state particle density by p^l r a (x). 
On the average, this 3D density will always be low in the parameter range 
considered here (in the sense that distance between particles is large com- 
pared to the 3D scattering length). The effective ID density can be either 
high or low, however. 

In parallel with the 3D Hamiltonian we consider the Hamiltonian for n 
Bosons in ID with delta interaction and coupling constant g > , i.e., 

H™ = jZ-d 2 /dz 2 ]+ g £ 5(z t - Zj ). (8.3) 
j=i i<*<i<^ 

We consider this Hamiltonian for the zj in an interval of length I in the 
thermodynamic limit, i — ► oo, n — > oo with p = njt fixed. The ground state 
energy per particle in this limit is independent of boundary conditions and 
can, according to |LL| . be written as 

el D (p)=p 2 e(g/p) , (8-4) 

with a function e(t) determined by a certain integral equation. Its asymp- 
totic form is e(t) ~ ht for t <C 1 and e(t) — > 7r 2 /3 for t — ► oo. Thus 

eh D (p)~l9P for g/p<l (8.5) 

and 

4 D (P) ~ (^/3)P 2 for <?/p>l ■ (8.6) 

This latter energy is the same as for non- interacting fermions in ID, which 
can be understood from the fact that (|8.3|) with g = oo is equivalent to a 
Hamiltonian describing free fermions. 
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Taking pe^f(p) as a local energy density for an inhomogeneous ID system 
we can form the energy functional 

/oo 
{\V^p(z)\ 2 + V L (z)p(z) + p(zfe(g/p{z))) dz . (8.7) 
-oo 

Its ground state energy is obtained by minimizing over all normalized den- 
sities, i.e., 



E m (N, L, g) =inf jffp] : p(z)>0, f 



p{z)dz = N 



Using convexity of the map p i— ► p 3 e(g/p), it is standard to show that there 
exists a unique minimizer of (jHUl) (see, e.g., jLSeYlj '). It will be denoted by 
PN,L,g- We also define the mean ID density of this minimizer to be 



1 

N 



(p N ,L,g(z)) 2 dz . (8.9) 



In a rigid box, i.e., for s = oo, p is simply N/L (except for boundary 
corrections), but in more general traps it depends also on g besides N and 
L. The order of magnitude of p in the various parameter regions will be 
described below. 

Our main result relates the 3D ground state energy of (|8.1|) to the ID 
density functional energy E 1T) (N,L,g) in the large TV" limit with g ~ a/r 2 
provided r/L and a/r are sufficiently small. To state this precisely, let e ± 
and b(x L ), respectively, denote the ground state energy and the normalized 
ground state wave function of — + V' L (x- L ). The corresponding quantities 
for —A 1 - + V r L (x ± ) are e^/r 2 and b r {x^) = (l/r)6(x _L /r). In the case that 
the trap is a cylinder with hard walls b is a Bessel function; for a quadratic 
V 1 - it is a Gaussian. 

Define g by 



87ra 



J \b[x L )\ A dx L = 8vra J \b r {x L )\ A dx L . (8.10) 



Our main result of this section is: 

Theorem 8.1 (From 3D to ID). Let N — > oo and simultaneously r/L — > 
and a/r — > in such a way that r 2 p ■ min{p, g] — » 0. Then 

E^ M (N,L,r,a) - Ne^/r 2 

An analogous result hold for the ground state density. Define the ID QM 
density by averaging over the transverse variables, i.e., 

-QM / \ f QM 



PlLJ^^)d^ ■ (8-12) 
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Let L := N/p denote the extension of the system in z-direction, and define 
the rescaled density p by 

Pa?l» = Jp(z/L) . (8.13) 

Note that, although p depends on N, L and g, \\p\\i = \\pW2 = 1, which 
shows in particular that L is the relevant scale in z-direction. The result for 
the ground state density is: 

Theorem 8.2 (ID limit for density). In the same limit as considered in 
Theorem \8.1i 

ft* (^pfi,rM) ~ = ( 8 - 14 ) 

in weak L 1 sense. 

Note that because of Q8.5|) and Q8.6|) the condition r 2 p ■ mm{p, g} — > is 
the same as 

e™(p) « 1/r 2 , (8.15) 

i.e., the average energy per particle associated with the longitudinal mo- 
tion should be much smaller than the energy gap between the ground and 
first excited state of the confining Hamiltonian in the transverse directions. 
Thus, the basic physics is highly quantum-mechanical and has no classical 
counterpart. The system can be described by a ID functional (|8.7|) . even 
though the transverse trap dimension is much larger than the range of the 
atomic forces. 

8.1. Discussion of the results. We will now give a discussion of the var- 
ious parameter regions that are included in the limit considered in Theo- 
rems 18.11 and 18.21 above. We begin by describing the division of the space 
of parameters into two basic regions. This decomposition will eventually be 
refined into five regions, but for the moment let us concentrate on the basic 
dichotomy. 

In SectionElwe proved that the 3D Gross-Pitaevskii formula for the energy 
is correct to leading order in situations in which N is large but a is small 
compared to the mean particle distance. This energy has two parts: The 
energy necessary to confine the particles in the trap, plus the internal energy 
of interaction, which is NAirap 3Y) . This formula was proved to be correct for 
a fixed confining potential in the limit N — » 00 with a 3 y0 3D — > 0. However, 
this limit does not hold uniformly if r/L gets small as N gets large. In other 
words, new physics can come into play as r/L — > and it turns out that 
this depends on the ratio of a/r 2 to the ID density, or, in other words, on 
g/p. There are two basic regimes to consider in highly elongated traps, i.e., 
when r <C L. They are 

• The ID limit of the 3D Gross-Pitaevskii regime 

• The 'true' ID regime. 
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The former is characterized by g/p <C 1, while in the latter regime g/p is 
of the order one or even tends to infinity. (If g/p — > oo the particles are 
effectively impenetrable; this is usually referred to as the Girardeau- Tonks 
region.) These two situations correspond to high ID density (weak interac- 
tion) and low ID density (strong interaction), respectively. Physically, the 
main difference is that in the strong interaction regime the motion of the 
particles in the longitudinal direction is highly correlated, while in the weak 
interaction regime it is not. Mathematically, this distinction also shows up 
in our proofs. The first region is correctly described by both the 3D and ID 
theories because the two give the same predictions there. That's why we 
call the second region the 'true' ID regime. 

In both regions the internal energy of the gas is small compared to the en- 
ergy of confinement. However, this in itself does not imply a specifically ID 
behavior. (If a is sufficiently small it is satisfied in a trap of any shape.) ID 
behavior, when it occurs, manifests itself by the fact that the transverse mo- 
tion of the atoms is uncorrelated while the longitudinal motion is correlated 
(very roughly speaking) in the same way as pearls on a necklace. Thus, the 
true criterion for ID behavior is that g/p is of order unity or larger and not 
merely the condition that the energy of confinement dominates the internal 
energy. 

We shall now briefly describe the finer division of these two regimes into 
five regions altogether. Three of them (Regions 1-3) belong to the weak 
interaction regime and two (Regions 4-5) to the strong interaction regime. 
They are characterized by the behavior of g/p as N —* oo. In each of these 
regions the general functional (|8.7j) can be replaced by a different, simpler 
functional, and the energy E 1P> (N, L, g) in Theorem 18. II bv the ground state 
energy of that functional. Analogously, the density in Theorem 18.21 can 
be replaced by the minimizer of the functional corresponding to the region 
considered. 

The five regions are 

• Region 1, the Ideal Gas case: In the trivial case where the interaction 
is so weak that it effectively vanishes in the large N limit and everything 
collapses to the ground state of —d 2 /dz 2 + V(z) with ground state energy 
ell, the energy E 1P> in Q8.11|) can be replaced by Ne^/L 2 . This is the case if 
g/p -C iV~ 2 , and the mean density is just p ~ N/L. Note that g/p <C N~ 2 
means that the 3D interaction energy per particle ~ ap 3D <C 1/L 2 . 

• Region 2, the ID GP case: In this region g/p ~ N~ 2 , with p ~ N/L. 
This case is described by a ID Gross-Pitaevskii energy functional of the form 

/oo 
{\V^-p{z)\ 2 + V L (z)p(z) + \gp(z) 2 ) dz , (8.16) 
-oo 

corresponding to the high density approximation (|8.5|) of the interaction 
energy in (|8.7j) . Its ground state energy, E GP , fulfills the scaling relation 
E GP (N, L, g) = NL- 2 E GP {1, l,NgL). 
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• Region 3, the ID TF case: iV~ 2 <C g/p -C 1, with p being of the 
order p ~ (N/L)(NgL)^ 1 ^ s+1 \ where s is the degree of homogeneity of the 
longitudinal confining potential V. This region is described by a Thomas- 
Fermi type functional 

/oo 
(V L (z)p(z) + \gp{zf) dz . (8.17) 
-oo 

It is a limiting case of Region 2 in the sense that NgL S> 1, but a/r 
is sufficiently small so that g/p -C 1, i.e., the high density approxima- 
tion in (|8.5|) is still valid. The explanation of the factor (NgL) 1 ^ 13 ^ 
is as follows: The linear extension L of the minimizing density of (|8.16|) 
is for large values of NgL determined by Vl{L) ~ g(N/L), which gives 
L ~ (N gL) l /( sJrl ^ L. In addition condition (|8.15|) requires gp -C r~ 2 , which 
means that Na/L(NgL) l / i ~ s+1 ^ « 1. The minimum energy of (|8.17[) has the 
scaling property E TF (N,L,g) = NL- 2 (NgL) s ^ s+1 ^ E TF (1, 1, 1). 

• Region 4, the LL case: g/p~l , with p ~ (iV/L)iV- 2 /( s+2 ), described 
by an energy functional 

/oo 
(F L (z)p(z)+p(z) 3 e( 5 /p(z)))cte. (8.18) 
-oo 

This region corresponds to the case g/p ~ 1, so that neither the high 
density (|8.5|) nor the low density approximation (|8.6j) is valid and the 
full LL energy 1)8.41) has to be used. The extension L of the system 
is now determined by Vi(L) ~ (N/L) 2 which leads to L ~ LN 2 ^ s+2 \ 
Condition 1)8.15(1 means in this region that Nr/L ~ N s ^ s+2 ^r/L — > 0. 
Since Nr/L ~ (p/ g)(a/r), this condition is automatically fulfilled if g/p 
is bounded away from zero and a/r — > 0. The ground state energy of (|8.18|) . 
E hL (N, L, g), is equal to Nj 2 E LL (l, 1, g/7), where we introduced the den- 
sity parameter 7 := (N / L)N- 2 ^ s+2 \ 

• Region 5, the GT case: g/p > 1, with p ~ (N/L)N- 2 ^ s+2 \ described 
by a functional with energy density ~ p 3 , corresponding to the Girardeau- 
Tonks limit of the LL energy density. It corresponds to impenetrable parti- 
cles, i.e, the limiting case g/p — > 00 and hence formula (|8.6j) for the energy 
density. As in Region 4, the mean density is here p ~ 7. The energy 
functional is 

/oo 
{V L (z)p(z) + (7r 2 /3)p(zf)dz , (8.19) 
-00 

with minimum energy i? GT (A r , L) = N^ 2 E GT (1, 1). 

As already mentioned above, Regions 1-3 can be reached as limiting cases 
of a 3D Gross-Pitaevskii theory. In this sense, the behavior in these regions 
contains remnants of the 3D theory, which also shows up in the fact that 
BEC prevails in Regions 1 and 2 (See |LSeY6j for details.) Heuristically, 
these traces of 3D can be understood from the fact that in Regions 1-3 the 
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ID formula for energy per particle, gp ~ aN/(r 2 L), gives the same result as 
the 3D formula, i.e., scattering length times 3D density. This is no longer 
so in Regions 4 and 5 and different methods are required. 

8.2. Outline of Proof. We now outline the main steps in the proof of 
Theorems 18 . 1 1 and 18 . 21 referring to |LSeY6j for full details. To prove Q8.11JI 
one has to establish upper and lower bounds, with controlled errors, on the 
QM many-body energy in terms of the energies obtained by minimizing the 
energy functionals appropriate for the various regions. The limit theorem 
for the densities can be derived from the energy estimates in a standard way 
by variation with respect to the external potential Vj_,. 

The different parameter regions have to be treated by different methods, 
a watershed lying between Regions 1-3 on the one hand and Regions 4-5 on 
the other. In Regions 1-3, similar methods as in the proof of the 3D Gross- 
Pitaevskii limit theorem discussed in Section El can be used. This 3D proof 
needs some modifications, however, because there the external potential was 
fixed and the estimates are not uniform in the ratio r/L. We will not go into 
the details here, but mainly focus on Regions 4 and 5, where new methods 
are needed. It turns out to be necessary to localize the particles by dividing 
the trap into finite 'boxes' (finite in z-direction), with a controllable particle 
number in each box. The particles are then distributed optimally among 
the boxes to minimize the energy, in a similar way as Eq. (|2.52[) was derived 
from Eq. (f2~47|) . 

A core lemma for Regions 4-5 is an estimate of the 3D ground state energy 
in a finite box in terms of the ID energy of the Hamiltonian (|8.3|) . I.e., we will 
consider the ground state energy of (|8.1j) with the external potential Vl{z) 
replaced by a finite box (in z-direction) with length £. Let E^ M (n, £, r, a) and 
E^ M (n, £, r, a) denote its ground state energy with Dirichlet and Neumann 
boundary conditions, respectively. 

Lemma 8.3. Let E^P(n,£,g) and E^f(n,£,g) denote the ground state en- 
ergy of \8JJ\) on L 2 ([0,£] n ), with Dirichlet and Neumann boundary condi- 
tions, respectively, and let g be given by \H.1()\) . Then there is a finite number 
C > such that 



E^ M (n,£,r,a) 




a\V8 



nr /a\ 

1 + 




Moreover, 



E^ M (n,£,r,a) 




(8.21) 



provided the term in square brackets is less than 1. 
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This Lemma is the key to the proof of Theorems 18.11 and 18.21 The reader 
interested in the details is referred to |LSeY6j . Here we only give a sketch 
of the proof of Lemma 18.31 

Proof of Lemma \H.'A We start with the upper bound (|8.21|) . Let ij) de- 
note the ground state of (|8.3[) with Dirichlet boundary conditions, normal- 
ly) 

ized by {ip\i/j) = 1, and let denote its two-particle density, normal- 
ized by f p^\z, z')dzdz' = 1. Let G and F be given by G{x\ , . . . , x n ) = 
ip(z 1 , . . .,z n )Uj =1 b r (xj) and F(xi, ...,x n ) = Yl^j f(Wi ~ Xj\)- Here / is 
a monotone increasing function, with < / < 1 and f(t) = 1 for t > R for 
some R > Ro. For t < R we shall choose f(t) = fo(t)/ fo(R), where /o is 
the solution to the zero-energy scattering equation for v a (|2.3j) . Note that 
fo(R) = 1 — a/R for R > Rq, and /g(i) < min{l, a/t}. We use as a trial 
wave function 

^(fi,...,x n ) = G(x 1 ,...,x n )F(x 1 ,...,aS n ) . (8.22) 

We first have to estimate the norm of ^. Using the fact that F is 1 
whenever no pair of particles is closer together than a distance R, we obtain 

, , , n(n — 1) irR 2 .,, ..a 

>1~ 9 - — b |. 8.23 

To evaluate the expectation value of the Hamiltonian, we use 

(tf | - Aj|*) = - J F 2 GAjG + J G 2 \VjF\ 2 (8.24) 
and the Schrodinger equation H n ^ g ip = E^ip. This gives 



MS 



<3 

2 



+ g 2 Ei^i'+E^-^D^ 



Now, since < / < 1 and /' > by assumption, F 2 < f(\xi 



x 




E|V i F| 2 <2^/'(|x,-f i |) 2 + 4 £ /'(Ixfe-fiD/'dxfc-Xjl) . (8.26) 

j=l i<j k<i<j 

Consider the first term on the right side of (|8.26j) . together with the last 
term in (|8.25j) . These terms are bounded above by 



n 



(n-1) J 6 r (f ± )\.(y ± )Vi 2) (^^)(/ / (l^-y1) 2 + ^a(|f-y|)/(|x-y|) 2 ). 

(8.27) 

Let 

h(z) = [ (/'(|f |) 2 + lv a (\x\)f(\x |) 2 ) ^ . (8.28) 
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Using Young's inequality for the integration over the _L-variables, we get 



n(n — 1) 



pf{z,z!)h{z- z')dzdz' 



\j: 2 



1.29) 



By similar methods, one can show that the contribution from the last term 
in l|8.26j) is bounded by 



3 

where 



n(n - l)(n - 2)- 



oo 



fe[|oo / P${z, z')k(z - z')dzdz , 



k(z) 



f'(\x\)d^ ■ 



3.30) 
$.31) 



Note that both h and k are supported in [— R, R]. 
Now, for any G H X (R), 



\\4>{z)\ 2 -\4>W\ 2 \<2\z-W( 



1/4 



d((> 



dz 



3/4 



3.32) 



(2) 

Applying this to p\ (z, z), considered as a function of z only, we get 



p^\z, z')h(z — z')dzdz — I h(z)dz I p [ ^ > (z,z)dz 



< 2R 1/2 / h(z)dz(^ 



dz\ 



\ 3/4 



(8.33) 



(2) 

where we used Schwarz's inequality, the normalization of p\ J and the sym- 
metry of ip. The same argument is used for 1)8.30(1 with h replaced by k. 

It remains to bound the second term in 1)8.25(1 . As in the estimate for the 
norm of ^, we use again the fact that F is equal to 1 as long as the particles 
are not within a distance R. We obtain 

n(n — 1) 



P% ,{z,z)dz I 1 



„2 H°ll4 



3.34) 



We also estimate g\n{n - 1) J p\'(z,z)dz < E^P and (ip\ - d 2 /dz 2 \ip) < 
E^P /n. We have f h(z)dz = 4-7ra(l — a/it!) -1 , and the terms containing k 
can be bounded by ||fe||oo < 27ra(l + ln(i?/a))/(l — a/r) and f k(z)dz < 
2itaR(l — a/(2R))/(l — a/r). Putting together all the bounds, and choosing 



R d 



ar 



(8.35) 



n 2 {l+gl) ' 
this proves the desired result. 

We are left with the lower bound (|8.2U|) . We write a general wave function 
* as 



n ) Y\ b r(s 



k I ) 



$.36) 



k=l 
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which can always be done, since b r is a strictly positive function. Partial 
integration gives 



i=i 



3, 3¥=i 



W K{xifdx k . 



k=l 



(8.37) 

Choose some R > Ro, fix i and Xj, j ^ i, and consider the Voronoi cell 
Clj around particle j, i.e., Qj = {x : \x — Xj\ < \x — x k \ for all k ^ j}. If 
Bj denotes the ball of radius R around Xj, we can estimate with the aid of 
Lemma 12.51 



b r (4) 2 (lVi/1 2 + \v a (\xi - Xj])]/] 2 ) dx 
min^gg brix- 1 -) 2 



> 



maxaf eS . 6 r (f J -) 2 j^ns. 



brixtyUQxi-XjDlfl 2 . (8.38) 



Here U is given in ()2.3U[) . For some 5 > let Bg be the subset of M 2 
where b(x ± ) 2 > S, and let xb s denote its characteristic function. Estimating 
max^gg. b r (x ) 2 < min^gg. b r (x^) 2 + 2(i?/r 3 )||V6 2 || 00 , we obtain 



min^gg. b r {x- 



u r-+±\2 - XBs@i /r) 1- 2- 
max^g B . b r (x^Y J V r 5 



R \\Vb 2 



1.39) 



Denoting k(i) the nearest neighbor to particle i, we conclude that, for < 
£< 1, 



jQ K(x^) 2 dx k 



k=l 



Vi/P + I^Uad^-^-DI/l 1 

3, 3¥=i 

>J2j [^if\ 2 + (l-e)\V l f\ 2 X mi n k \z t -z k \> R (z l ) 

i=l 

n 

+a'U(\xi - x k{i) \) X B s (4(i)M\f\ 2 } II b r(4) 2 dx k , (8.40) 



k=l 



where a' = a(l - e)(1 - 2i2||V6 2 || 00 /rJ). 
Define i^i, . . . , z n ) > by 



/n 
k=i 

Neglecting the kinetic energy in _L-direction in the second term in (J8.4(J|) 
and using the Schwarz inequality to bound the longitudinal kinetic energy 
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of / by the one of F, we get the estimate 

(m\*)- — > 

n „ n 

-z k \>R{zi) \ \\dz k 
1=1 J k=l 

n „ n 

+e / [ £ i v ^/i 2 + a Mi^-%)Dx%(4 w / r )i/i 2 ]n^(^) 2d ^' 

i=l ^ k=l 

(8.42) 

where <9j = d/dzj, and V -1- denotes the gradient in _L-direction. We now 
investigate the last term in l|8.42|) . Consider, for fixed Zi,...,z n , the expres- 
sion 

n „ n 

E / H^/l' + a'C/d^-^DxB.^w/OI/l 2 ] Hb r (4) 2 d4 . 

i=l fe=l 

(8.43) 

To estimate this term from below, we use Temple's inequality, as in 
Subsect. 12.21 Let e ± denote the gap above zero in the spectrum of 
— A 1 - + V ± — e ± , i.e., the lowest non-zero eigenvalue. By scaling, e^/r 2 is 
the gap in the spectrum of —A 1 - + — e ± /r 2 . Note that under the trans- 
formation 4> l—> b~ l (p this latter operator is unitarily equivalent to V^* • V -1- 
as an operator on L 2 (M 2 , b^x^) 2 dx^) , as considered in ()8.43j) . Hence also 
this operator has e^/r 2 as its energy gap. Denoting 

// n \ k n 

[J2 U (\^- S k(i)\)XB s (4({)M) l\b r (4) 2 dxi , (8.44) 
\i=l / k=l 

Temple's inequality implies 

633 > \f W(U) (i - °-^ ^ /r2 1 _„, w ) ■ ("5) 

Now, using (l2~3uTl and Schwarz's inequality, {U 2 ) < 3n(i? 3 - J2g)— 1 <I7>, and 
Therefore 

|2„// 



I1SJ) > |.F|V(L0 , (8.47) 

where we put all the error terms into the modified coupling constant a" . It 
remains to derive a lower bound on ([/). Let 

d(z-z')= [ b r {x^fb r ^) 2 U{\x-y\)xB s {t/r)dx^dt ■ (8.48) 
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Note that d(z) = if \z\ > R. An estimate similar to (|2.3(i|) gives 

(U)>^d(z t -z j )(l-(n-2)^\\b\\l Q ) . (8.49) 

Note that, for an appropriate choice of R, d is close to a ^-function with the 
desired coefficient. To make the connection with the 5- function, we can use 
a bit of the kinetic energy saved in (|8.42j) to obtain 



n — 1 



\diF\ z + a'"d(zi- Zj )\F\ 



dzi 

1/2^ 



> y max \F\\ [R/ _ R] { Zj ) ( 1 - (^LAj/r) ) . (8.50) 
\zi-zj\<R v \ e j j 

Putting all the previous estimates together, we arrive at 

<*|2r|¥> - > ^ / [(i-£)|^| 2 Xminfek _, fel >^]n^ 

' i=l J k=l 

+ Y^y / , ma xj F \ 2 xw-R]( z j) n ^ 



(8.51) 

for an appropriate coupling constant g" that contains all the error terms. 
Now assume that (n + 1)R < t. Given an F with J \F\ 2 dzi • ■ ■ dz n = 1, 
define, for < z\ < Z2 < ■ ■ ■ < z n < £ — (n + l)R, 

i>(z u . ..,z n ) = F( Zl + R,z 2 + 2R, z 3 + 3R,...,z n + nR) , (8.52) 

and extend the function to all of [0, 1 — (n+ l)R] n by symmetry. A simple 
calculation shows that 

JEHD > {il>\H'\il>) > (l-e)E^(n,e-(n + l)R,g")(i>\i>) 

> (l-£)4 D (n,l )5 ")(#) , (8.53) 

where H' is the Hamiltonian (|8.3[) with a factor (1 — e) in front of the kinetic 
energy term. 

It remains to estimate (ip\ip). Using that F is related to the true ground 
state \I/ by (|8.41f) , we can estimate it in terms of the total QM energy, namely 



<V#> > l-^(^ M (iV,^,r,a) 



ne~ 



1/2 



-2n^-4nR(-E^ M (n,e,r,a)-^\ . (8.54) 

Collecting all the error terms and choosing 

/a\i/4 /a\i/8 p /a\i/8 

R = r (r) ' £= (r) ' (8 - 55) 

(|8.53f) and (|8.54l) lead to the desired lower bound. ■ 
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As already noted above, Lemma 18.31 is the key to the proof of Theo- 
rems 18.11 and 18.21 The estimates are used in each box, and the particles 
are distributed optimally among the boxes. For the global lower bound, 
superadditivity of the energy and convexity of the energy density p 3 e(g/p) 
are used, generalizing corresponding arguments in Section |2j We refer to 
jhSeYfij for details. 

9. The Charged Bose Gas, the One- and Two-Component Cases 

The setting now changes abruptly. Instead of particles interacting with a 
short-range potential v(\x{ — Xj\) they interact via the Coulomb potential 



v(\Xi - Xj\) = \Xi 



1-1 



(in 3 dimensions). The unit of electric charge is 1 in our units. 

We will here consider both the one- and two-component gases. In the 
one-component gas (also referred to as the one-component plasma or bosonic 
jellium) we consider positively charged particles confined to a box with a 
uniformly charged background. In the two-component gas we have particles 
of both positive and negative charges moving in all of space. 

9.1. The One-Component Gas. In the one-component gas there are N 
positively charged particles in a large box A of volume L 3 as before, with 
p = N/L 3 . 

To offset the huge Coulomb repulsion (which would drive the particles to 
the walls of the box) we add a uniform negative background of precisely the 
same charge, namely density p. Our Hamiltonian is thus 

N 

H$ =^-A*Ai-y(x f )+ J2 <Wi-Sj\) + C (9.1) 

i=l l<i<j<N 

with 

V(x) = P I \x — y\~ 1 dy and C = -p f V(x)dx . 

J A 2 JA 

We shall use Dirichlet boundary conditions. As before the Hamiltonian acts 
on symmetric wave functions in L 2 (A N , dx% ■ ■ ■ dxjy). 

Each particle interacts only with others and not with itself. Thus, despite 
the fact that the Coulomb potential is positive definite, the ground state 
energy E$ can be (and is) negative (just take =const.). This time, large 
p is the 'weakly interacting' regime. 

We know from the work in |LN| that the thermodynamic limit eo(p) de- 
fined as in 1)2. 2j) exists. It also follows from this work that we would, in 
fact, get the same thermodynamic energy if we did not restrict the number 
of particles N, but considered the grand-canonical case where we minimize 
the energy over all possible particle numbers, but keeping the background 
charge p fixed. 
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Another way in which this problem is different from the previous one is 
that perturbation theory is correct to leading order. If one computes (\l/, H^) 
with \I> =const, one gets the right first order answer, namely 0. It is the next 
order in 1/p that is interesting, and this is entirely due to correlations. In 
1961 Foldy [F] calculated this correlation energy according to the prescrip- 
tion of Bogolubov's 1947 theory. That theory was not exact for the dilute 
Bose gas, as we have seen, even to first order. We are now looking at second 
order, which should be even worse. Nevertheless, there was good physical 
intuition that this calculation should be asymptotically exact. Indeed it is, 
as proved in |LSoj and |Soj . 

The Bogolubov theory states that the main contribution to the energy 
comes from pairing of particles into momenta k,—k and is the bosonic ana- 
logue of the BCS theory of superconductivity which came a decade later. 
I.e., is a sum of products of terms of the form expji/c • (Si — xj)}. 

The following theorem is the main result for the one-component gas. 

Theorem 9.1 (Foldy's law for the one-component gas). 



This is the first example (in more than 1 dimension) in which Bogolubov's 
pairing theory has been rigorously validated. It has to be emphasized, how- 
ever, that Foldy and Bogolubov rely on the existence of Bose-Einstein con- 
densation. We neither make such a hypothesis nor does our result for the 
energy imply the existence of such condensation. As we said earlier, it is 
sufficient to prove condensation in small boxes of fixed size. 

Incidentally, the one-dimensional example for which Bogolubov's theory 
is asymptotically exact to the first two orders (high density) is the repulsive 
delta-function Bose gas |LL| . for which there is no Bose-Einstein condensa- 
tion. 

To appreciate the —p l / A nature of (|9.2j) . it is useful to compare it with 
what one would get if the bosons had infinite mass, i.e., the first term in 
()9.1j) is dropped. Then the energy would be proportional to —p 1 ^ 3 as shown 
in |LN| . Thus, the effect of quantum mechanics is to lower 5 to 7. 

A problem somewhat related to bosonic jellium is fermionic jellium. Graf 
and Solovej |GSj have proved that the first two terms are what one would 
expect, namely 



where Ctf is the usual Thomas-Fermi constant and Cd is the usual Dirac 
exchange constant. 

It is supposedly true, for both bosonic and fermionic particles, that there 
is a critical mass above which the ground state should show crystalline or- 
dering (Wigner crystal), but this has never been proved and it remains an 
intriguing open problem, even for the infinite mass case. A simple scaling 




(9.2) 



e (p) = C TFP 5 / 3 -C D p i / 3 + o(p 4 / 3 ) 



(9.3) 
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shows that large mass is the same as small p, and is thus outside the region 
where a Bogolubov approximation can be expected to hold. 

As for the dilute Bose gas, there are several relevant length scales in 
the problem of the charged Bose gas. For the dilute gas there were three 
scales. This time there are just two. Because of the long range nature of the 
Coulomb problem there is no scale corresponding to the scattering length a. 
One relevant length scale is again the interparticle distance p~ 1 ^ . The other 
is the correlation length scale £ cor ~ p -1 / 4 (ignoring the dependence on p). 
The order of the correlation length scale can be understood heuristically as 
follows. Localizing on a scale £ cor requires kinetic energy of the order of £~ 2 r . 
The Coulomb potential from the particles and background on the scale £ cor 
is (p£^ OT )/£ COT . Thus the kinetic energy and the Coulomb energy balance 
when £ cor ~ p~ l / A . This heuristics is however much too simplified and hides 
the true complexity of the situation. 

Note that in the high density limit £ cor is long compared to the interpar- 
ticle distance. This is analogous to the dilute gas where the scale £ c is also 
long compared to the interparticle distance [see (|2,12|) ]. There is however 
no real analogy between the scale £ cov for the charged gas and the scale £ c 
for the dilute gas. In particular, whereas eo(p) for the dilute gas is, up to 
a constant, of the same order as the kinetic energy ~ pt^ 1 we have for the 
charged gas that eo(p) </> £~ 2 r = p 1 ! 2 . The reason for this difference is that 
on average only a small fraction of the particles in the charged gas actually 
correlate. 

9.2. The Two- Component Gas. Now we consider N particles with char- 
ges ±1. The Hamiltonian is thus 

N 



This time we are interested in Eq (N) the ground state energy of H N 
minimized over all possible combination of charges = ±1, i.e., we do not 
necessarily assume that the minimum occurs for the neutral case. Restricting 
to the neutral case would however not change the result we give below. 

(2) 

An equivalent formulation is to say that Eq (N) is the ground state 
energy of the Hamiltonian acting on all wave functions of space and charge, 
i.e., functions in I? ((M 3 x {—1,1})^). As mentioned in the introduction, 
and explained in the beginning of the proof of Thm. 12.21 for the calculation 
of the ground state energy we may as usual restrict to symmetric functions 
in this Hilbert space. 

For the two-component gas there is no thermodynamic limit. In fact, 
Dyson [D2j proved that E^\n) was at least as negative as - (const) N 7 / 5 
as N — > oo. Thus, thermodynamic stability (i.e., a linear lower bound) 
fails for this gas. Years later, a lower bound of this —N 7 / 5 form was finally 
established in |CLY| . thereby proving that this law is correct. 



H' 
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The connection of this —N 7 ^ law with the jellium — p 1//4 law (for which 
a corresponding lower bound was also given in |CLY| ) was pointed out by 
Dyson |D2j in the following way. Assuming the correctness of the — p 1 / 4 
law, one can treat the 2-component gas by treating each component as a 
background for the other. What should the density be? If the gas has 
a radius L and if it has N bosons then p = NL~ 5 . However, the extra 
kinetic energy needed to compress the gas to this radius is NL~ 2 . The total 
energy is then NL~ 2 — N p 1//4 , and minimizing this with respect to L gives 
L ~ iV -1 / 5 and leads to the — iV 7 / 5 law. The correlation length scale is now 

In |D2j Dyson conjectured an exact asymptotic expression for (N) 



for large N. That this asymptotics, as formulated in the next theorem, is 
indeed correct is proved in |LSo2j and jSoj . 

Theorem 9.2 (Dyson's law for the two-component gas). 

J im \flP-= M \f" I |V*Mo/ * 5/2 0<d>, [ <& 2 = l), (9.4) 

Ar_oo jV 7 / 5 [ J R 3 J R 3 J R 3 J 

where Iq is the constant from Foldy 's law: 

_ 2 r(3/4) 

°"5f(574)W ' 

This asymptotics can be understood as a mean field theory for the gas 
density, very much like the Gross-Pitaevskii functional for dilute trapped 
gases, where the local energy described by Foldy's law should be balanced 
by the kinetic energy of the gas density. Thus if we let the gas density be 
given by 4> 2 then the "mean field" energy should be 

|V0| 2 -/ O f 5/2 - (9.5) 



Here fcf) 2 = N. If we now define $(x) = iV~ 8 / 5 (/)(A r ~ 1//5 f ) we see that 
J <J> 2 = 1 and that the above energy is 

N 7 / 5 (fM [ |V$| 2 -hi <K 5 / 2 

It may be somewhat surprising that it is exactly the same constant Iq 
that appears in both the one- and two-component cases. The reason that 
there are no extra factors to account for the difference between one and two 
components is, as we shall see below, a simple consequence of Bogolubov's 
method. The origin of this equivalence, while clear mathematically, does 
not appear to have a simple physical interpretation. 

9.3. The Bogolubov Approximation. In this section we shall briefly ex- 
plain the Bogolubov approximation and how it is applied in the case of the 
charged Bose gas. The Bogolubov method relies on the exact diagonalization 
of a Hamiltonian, which is quadratic in creation and annihilation operators. 



THE QUANTUM-MECHANICAL MANY-BODY PROBLEM: THE BOSE GAS 69 

For the charged Bose gas one only needs a very simple case of the general 
diagonalization procedure. On the other hand, the operators that appear 
are not exact creation and annihilation operators. A slightly more general 
formulation is needed. 

Theorem 9.3 (Simple case of Bogolubov's method). 

Assume that b± ± are four (possibly unbounded) commuting operators satis- 
fying the operator inequality 

[h,eiK,e] ^ 1 f or al1 e,r = ±. 
Then for all real numbers A,B + ,B- > we have 

r,e=±l 

+ ]T ^KB~iee'(b* +<e b +:e , + 6*_ ie 6_ je , + b* + , e b*_, e , + & + , e &_, e ,) 

e,e'=±l 

> -(A + B+ + B-) + y/(A + B+ + BJf - {B+ + BJf. 
If 6 ± ± are four annihilation operators then the lower bound is sharp. 
Proof. Let us introduce 

4 = (B + +B-)- 1/2 (B 1 / 2 b* ±;+ - B]/ 2 b*±_), 

and 

4 = (B+ + B-)- l/2 {B l _! 2 b* ±)+ +B 1 l 2 b* ±) _). 
Then these operators satisfy 

[d + ,d* + ]<l, [d_,d*_]<l. 
The operator that we want to estimate from below may be rewritten as 

A{d* + d + + d*_d_+ c* + c + +c*_c_) 

+ (B + + B-) (d* + d + + d*_d_ + d* + d*_ + d + d_) . 

We may now complete the squares to write this as 

A(c* +C+ + c*_c_) + D(d* + + Xd_)(d* + + Xd_)* 

+D(d*_ + Xd + )(d*_ + Xd + )* - DX 2 ([d + ,d* + ] + [d_,d*_]) 

if 

D(l + A 2 ) =A + B+ +B-, 2DX = B+ + B-. 
We choose the solution A = 1 + B+ + B ^ (l + j j — 1. Hence 

DX 2 = \(a + B+ + B-- V(A + B+ + B-) 2 - (B+ + £_) 2 ) . 
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In the theorem above one could of course also have included linear terms in 
6 Tj e in the Hamiltonian. In the technical proofs in iLSo1 lLSo2] the Bogolubov 
diagonalization with linear terms is indeed being used to control certain error 
terms. Here we shall not discuss the technical details of the proofs. We have 
therefore stated the theorem in the simplest form in which we shall need it 
to derive the leading contribution. 

In our applications to the charged Bose gas the operators b± e will cor- 
respond to the annihilation of particles with charge e = db and momenta 
±fc for some k E M 3 . Thus, only equal and opposite momenta couple. In a 
translation invariant case this would be a simple consequence of momentum 
conservation. The one-component gas is not translation invariant, in our 
formulation. The two-component gas is translation invariant, but it is nat- 
ural to break translation invariance by going into the center of mass frame. 
In both cases it is only in some approximate sense that equal and opposite 
momenta couple. 

In the case of the one-component gas we only need particles of one sign. 
In this case we use the above theorem with b±- = and B- = 0. 

We note that the lower bounds in Theorem 19.31 for the one- and two- 
component gases are the same except for the replacement of B+ in the one- 
component case by £>+ + B- in the two-component case. In the application 
to the two-component gas i3+ and £>_ will be proportional to the particle 
densities for respectively the positive or negatively charged particles. For 
the one-component gas B + is proportional to the background density. 

The Bogolubov diagonalization method cannot be immediately applied 
to the operators H]p or since these operators are not quadratic in 

creation and annihilation operators. In fact, they are quartic. They have 
the general form 

^tapa^ap + \ ^2 w a/3 ^a* a a}a u a^ (9.6) 

with 

t a p = (a\T\(3), w a p^ = (a(3\W\fj,v), 

where T is the one-body part of the Hamiltonian and W is the two-body-part 
of the Hamiltonian. 

The main step in Bogolubov's approximation is now to assume Bose- 
Einstein condensation, i.e., that almost all particles are in the same one- 
particle state. In case of the two-component gas this means that almost 
half the particles are positively charged and in the same one-particle state 
as almost all the other half of negatively charged particles. We denote this 
condensate state by the index a = in the sums above. Based on the 
assumption of condensation Bogolubov now argues that one may ignore all 
terms in the quartic Hamiltonian above which contain 3 or 4 non-zero indices 
and at the same time replace all creation and annihilation operators of the 
condensate by their expectation values. The result is a quadratic Hamilton- 
ian (including linear terms) in the creation and annihilation with non-zero 
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index. This Hamiltonian is of course not particle number preserving, reflect- 
ing the simple fact that particles may be created out of the condensate or 
annihilated into the condensate. 

In Section EH below it is explained how to construct trial wave functions 
for the one- and two-component charged gases whose expectations agree 
essentially with the prescription in the Bogolubov approximation. The de- 
tails are to appear in jSoj . This will imply upper bounds on the energies 
corresponding to the asymptotic forms given in Theorems 19.11 and 19 . 21 

In |LSo| ILSo2j it is proved how to make the steps in the Bogolubov ap- 
proximation rigorous as lower bounds. The main difficulty is to control the 
degree of condensation. As already explained it is not necessary to prove 
condensation in the strong sense described above. We shall only prove con- 
densation in small boxes. Put differently, we shall not conclude that most 
particles are in the same one-particle state, but rather prove that most par- 
ticles occupy one-particle states that look the same on short scales, i.e., that 
vary slowly. Here the short scale is the correlation length scale £ cor . 

9.4. The Rigorous Lower Bounds. As already mentioned we must local- 
ize into small boxes of some fixed size £. This time we must require £ cor <C £■ 
For the one-component gas this choice is made only in order to control the 
degree of condensation. For the two-component gas it is required both to 
control the order of condensation, and also to make a local constant density 
approximation. The reason we can control the degree of condensation in a 
small box is that the localized kinetic energy has a gap above the lowest 
energy state. In fact, the gap is of order £~ 2 . Thus if we require that £ is 
such that Nl~ 2 is much greater than the energy we may conclude that most 
particles are in the lowest eigenvalue state for the localized kinetic energy. 
We shall always choose the localized kinetic energy in such a way that the 
lowest eigenstate, and hence the condensate, is simply a constant function. 

9.4.1. Localizing the interaction. In contrast to the dilute gas the long range 
Coulomb potential prevents us from simply ignoring the interaction between 
the small boxes. To overcome this problem we use a sliding technique first 
introduced in |( "LY| . 

Theorem 9.4 (Controlling interactions by sliding). Let X be a smooth 
approximation to the characteristic function of the unit cube centered at the 
origin. Fori > and z G M 3 letXz(x) = x{{x — z)/£). There exists anu > 
depending on X (in such a way that it tends to infinity as X approximates 
the characteristic function) such that 




f* 2 Yf E 




Nuj 
~2l 



for all x\, . . . £l 3 and e±, . . . 

wg(x, y) 



Xs{x)Y u /i{x - y)Xz(y) 



±1, where 
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with Y^{x) = \x\ x exp(— p\x\) being the Yukawa potential. 

The significance of this result is that the two-body potential wg is localized 
to the cube of size £ centered at £z. The lower bound above is thus an integral 
over localized interactions sliding around with the integration parameter. 

We have stated the sliding estimate in the form relevant to the two- 
component problem. There is an equivalent version for the one-component 
gas, where the sum of the particle-particle, particle-background, and back- 
ground-background interactions may be bounded below by corresponding 
localized interactions. 

Since I 3> £ cor the error in the sliding estimate is much smaller than 
toN/£ cor , which for both the one and two-component gases is of order u> 
times the order of the energy. Thus, since £ is much bigger than £ COT , we 
have room to let w be very large, i.e., X is close to the characteristic function. 

9.4.2. Localizing the kinetic energy. Having described the technique to con- 
trol the interaction between localized regions we turn next to the localization 
of the kinetic energy. 

For the two-component gas this is done in two steps. As already men- 
tioned it is natural to break the translation invariance of the two-component 
gas. We do this by localizing the system into a box of size L' S> iV -1 / 5 (which 
as we saw is the expected size of the gas) as follows. By a partition of unity 
we can divide space into boxes of this size paying a localization error due 
to the kinetic energy of order NL'~ 2 <C iV 7 / 5 . We control the interaction 
between these boxes using the sliding technique. 

We may now argue, as follows, that the energy is smallest if all the par- 
ticles are in just one box. For simplicity we give this argument for the case 
of two boxes. Suppose the two boxes have respective wave functions ip and 
ip. The total energy of these two non-interacting boxes is E + E. Now 
put all the particles in one box with the trial function ^> = tpip. The fact 
that this function is not bosonic (i.e., it is not symmetric with respect to all 
the variables) is irrelevant because the true bosonic ground state energy is 
never greater than that of any trial state (Perron-Frobenius Theorem) . The 
energy of Vl/ is 



where p^ and p^ are the respective charge densities of the states ip and ip. 
We claim that the last Coulomb term can be made non-positive. How? If 
it is positive then we simply change the state ip by interchanging positive 
and negative charges (only in ip and not in ip). The reader is reminded that 
we have not constrained the number of positive and negative particles but 
only their sum. This change in ip reverses the relative charge of the states 
ip and ip so, by symmetry the energies E and E do not change, whereas the 
Coulomb interaction changes sign. 



E + E + 
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The localization into smaller cubes of size £ can however not be done by a 
crude partition of unity localization. Indeed, this would cost a localization 
error of order N£~ 2 , which as explained is required to be of much greater 
order than the energy. 

For the one-component charged gas we may instead use a Neumann lo- 
calization of the kinetic energy, as for the dilute Bose gas. If we denote by 

(z) -> 

A) ' the Neumann Laplacian for the cube of size i centered at z we may, in 
the spirit of the sliding estimate, write the Neumann localization Laplacian 
in all of M. 3 as 

" -AfW. 

In order to write the localized kinetic energy in the same form as the localized 
interaction we must introduce the smooth localization X as in Theorem 19.41 
This can be achieved by ignoring the low momentum part of the kinetic 
energy. 

More precisely, there exist e{x) and s(x) such that e(x) — * and s(x) — ► 
as X approaches the characteristic function of the unit cube and such that 
(see Lemma 6.1 in jLSop 

-Af > (1 - e(x))VzXz(£)F Hx) (-A)Xg{$)r? (9.7) 

where denotes the projection orthogonal to constants in the cube of size 
i centered at z and 



F s (u) 



u 2 



u + s 2 



For u <C s~ 2 we have that F s (u) <C u. Hence the effect of F in the operator 
estimate above is to ignore the low momentum part of the Laplacian. 

For the two-component gas one cannot use the Neumann localization 
as for the one-component gas. Using a Neumann localization ignores the 
kinetic energy corresponding to long range variations in the wave function 
and one would not get the kinetic energy term J /j,\X7Q\ 2 in Q9.4|) . This 
is the essential difference between the one- and two-component cases. This 
problem is solved in |LSo2j where a new kinetic energy localization technique 
is developed. The idea is again to separate the high and low momentum part 
of the kinetic energy. The high momentum part is then localized as before, 
whereas the low momentum part is used to connect the localized regions by 
a term corresponding to a discrete Laplacian. (For details and the proof the 
reader is referred to |LSo2j .) 

Theorem 9.5 (A many body kinetic energy localization). Let Xz> 

Vg and F s be as above. There exist e{x) and s(x) such that e{x) — ► 
and s(x) — ► as X approaches the characteristic function of the unit cube 
and such that for all normalized symmetric wave functions \P in L 2 
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{-1, 1})^) and all UcR 3 we have 



N 



i=l 





/[ 




Jn 1 



Y^(S*(t{z + y))-S*(lz)f 



+±r 2 

yei? 
Ial=i 

-const. r 2 Vol(^) 



dz 



where 



S*(z) = J (a* 0+ (z)a 0+ (z) + a5_(z)a _ + 1 



w;zf/j a 0= |-(2) foem^ i/ie annihilation of a particle of charge ± in the state given 
by the normalized characteristic function of the cube of size £ centered at z. 

The first term in the kinetic energy localization in this theorem is the 
same as in (|9.7j) . The second term gives rise to a discrete Laplacian for the 
function S^(£z), which is essentially the number of condensate particles in 
the cube of size £ centered at Iz. Since we will eventually conclude that 
most particles are in the condensate this term will after approximating the 
discrete Laplacian by the continuum Laplacian lead to the term J /i|V</>| 2 in 
Q9.5JI . We shall not discuss this any further here. 

When we apply this theorem to the two-component gas the set £0, will 
be the box of size U discussed above. Hence the error term £~ 2 Vol(£l) will 
be of order L l3 /£~ 5 <C (N' 2 / 5 £)- 5 (N 1 / B L , ) 3 N 7 / 5 . Thus since £ > N^ 5 we 
may still choose L' 3> iV^ 1 / 5 , as required, and have this error term be lower 
order than iV 7 / 5 . 

9.4.3. Controlling the degree of condensation. After now having localized 
the problem into smaller cubes we are ready to control the degree of con- 
densation. We recall that the condensate state is the constant function in 
each cube. Let us denote by rig the number of excited (i.e., non-condensed 
particles) in the box of size £ centered at z. Thus for the two-component gas 
aQ + (z)a Q+ (z) + aQ_(z)a _(z) is the total number of particles in the box 
and a similar expression gives the particle number for the one-component 
gas. 

As discussed above we can use the fact that the kinetic energy localized 
to a small box has a gap above its lowest eigenvalue to control the number 
of excited particles. Actually, this will show that the expectation (^,n^) 
is much smaller than the total number of particles in the box for any state 
^ with negative energy expectation. 

One needs, however, also a good bound on (^,^1^) to control the 
Coulomb interaction of the non-condensed particles. This is more difficult. 
In |LSoj this is not achieved directly through a bound on (^,ng^) in the 
ground state. Rather it is proved that one may change the ground state 
without changing its energy very much, so that it only contains values of rig 
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localized close to (^,ng^/). The following theorem gives this very general 
localization technique. Its proof can be found in jLSo . 

Theorem 9.6 (Localizing large matrices). Suppose that A is an N + 
1 x N + 1 Hermitean matrix and let A k , with k = 0,1,..., N, denote the 
matrix consisting of the k th supra- and infra- diagonal of A. Let tj) G C N+1 
be a normalized vector and set d}. = (tl),A k tj) and A = (ip,Aip) = X^fcLo4- 
(tj) need not be an eigenvector of A.) 

Choose some positive integer M < N + 1. Then, with M fixed, there is 
some n G [0, N + 1 — M] and some normalized vector (j> G C N+1 with the 
property that (j)j = unless n + \<j<n + M (i.e., 4> has length M) and 
such that 

„ M-l N 

(4>,A^<X + W Y1 k 2 \d k \+CY, 141 , (9-8) 

k=l k=M 

where C > is a universal constant. (Note that the first sum starts with 
k = I.) 

To use this theorem we start with a ground state (or approximate ground 
state) $ to the many body problem. We then consider the projections of 
onto the eigenspaces of rig. Since the possible eigenvalues run from to N 
these projections span an at most N + 1 dimensional space. 

We use the above theorem with A being the many body Hamiltonian 
restricted to this N + 1 dimensional subspace. Since the Hamiltonian can 
change the number of excited particles by at most two we see that dk van- 
ishes for k > 3. We shall not here discuss the estimates on d\ and c?2 (see 
|LSo| ILSo2j ) . The conclusion is that we may, without changing the energy 
expectation of ^ too much, assume that the values of ng run in an interval 
of length much smaller than the total number of particles. We would like 
to conclude that this interval is close to zero. This follows from the fact 
that any wave function with energy expectation close to the minimum must 
have an expected number of excited particles much smaller than the total 
number of particles. 

9.4.4. The quadratic Hamiltonian. Using our control on the degree of con- 
densation it is now possible to estimate all unwanted terms in the Hamil- 
tonian, i.e., terms that contain 3 or more creation or annihilation operators 
corresponding to excited (non-condensate) states. The proof which is a 
rather complicated bootstrapping argument is more or less the same for the 
one- and two-component gases. The result, in fact, shows that we can ig- 
nore other terms too. In fact if we go back to the general form (|9.6j) of the 
Hamiltonian it turns out that we can control all quartic terms except the 
ones with the coefficients: 



U>a/300, U)00a/3, and Wo a /30- 
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To be more precise, let u a , a = 1, . . . be an orthonormal basis of real func- 
tions for the subspace of functions on the cube of size i centered at z or- 
thogonal to constants, i.e, with vanishing average in the cube. We shall now 
omit the subscript z and let ao± be the annihilation of a particle of charge 
±1 in the normalized constant function in the cube (i.e., in the condensate). 
Let a a ± with q / be the annihilation operator for a particle of charge 
±1 in the state u a . We can then show that the main contribution to the 
localized energy of the two-component gas comes from the Hamiltonian 

oo 

#local = ^ taf3a*ae a f3e 

a, 0=1 
e=±l 

+ \ ^ ee ' W al3( 2a 0e a *ae> a 0e' a f3e + a *0e a *0e' a a e' a /3e + a *ae a }e' a 0e' a 0e), 

a, 0=1 
e,e'=±l 



where 



and 



t a p = fi(u a ,VgXz(x)F es ( x) (-A)Xz(x)'Pgup) 



a/ 3=£ 3 J j u a (x)Xz{x)Y u / e (x - y)Xz{y)uf3{x)dxdy. 



w, 

In i?iocai we have ignored all error terms and hence also e(x) ~ and 

/x 2 «i. 

In the case of the one-component gas we get exactly the same local Hamil- 
tonian, except that we have only one type of particles, i.e, we may set 
a a _ = above. 

Let v± = Yla=o a *a± a a± be the total number of particles in the box with 
charge ±1. For k € M? we let X^g{x) = Xs{x)e lk ' x . We then introduce the 
operators 

b %± = (i 3 ^±r 1/2 a ± (V,X l :Ja* 0± , 

where a±(V?Xj: -) = SaLiC^n-i ««)a a ± annihilates a particle in the state 
X^ g with charge ±1. It is then clear that the operators b^ ± all commute 
and a straightforward calculation shows that 

h±> b W ^ (^ 3 ^)- 1 ||^|| 2 ao±«o± < 1- 
If we observe that 

oo 



t aP a* ae a pe = (2vr)- 3 / ^F la{x) {P)^a e {V 2 X^)*a e {V s X^ s )dk 

*,0=i J e=± 

5=±1 

> (2 7 r)-^/^(x)(^)E^e' 
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we see that 



J e=± 

+ £ t/MV^ee'ibtb^, + b*_ % b_^, + bl/_ %el + b_ Ae ,)dk 

ee'=± 
a/3=l 

The last term comes from commuting Oo± a o± to a o± a o±- ^ ^ s eas y to see 
that this last term is a bounded operator with norm bounded by 

const. {v + + )€~ 3 1| 1^ /£ ||oo < const. uj~ 2 (v + + vJ)t~ x . 

When summing over all boxes we see that the last term above gives a con- 
tribution bounded by const. io~ 2 Nt~ x = uj^ 2 (N 2 ^IY 1 N 7 ^ which is lower 
order than the energy. 

The integrand in the lower bound on H\ oc& \ is precisely an operator of the 
form treated in the Bogolubov method Theorem 19.31 Thus up to negligible 
errors we see that the operator H\ OC3 x is bounded below by 



2 

where 



i(2vr)- 3 J -(A(k) + B(k)) + y/(A(k) + B(k)) 2 -B(k) 2 dk, 



A(k) = ^e 3 F £sM (k 2 ) and B(k) = vY w/e (k) 



with v = v + + V- being the total number of particles in the small box. A 
fairly simple analysis of the above integral shows that we may to leading 
order replace A by /i£ 3 k 2 and B(k) by ^v\k\~ 2 , i.e., we may ignore the cut- 
offs. The final conclusion is that the local energy is given to leading order 
by 



— ^ f ^v\k\~ 2 + n£ 3 \k\ 2 - J (4vri/|£|- 2 + Lil 3 \k\ 2 ) 2 - (4vr^|fc|- 2 ) 2 dk 

2(277)- J V 

= _2 1 / 2 vr -3/4 zy (JL\ V4 j'^l + x 4 - x 2 (2 + x 4 ) 1 / 2 dx. 



1/4 

^ f — 1 

If we finally use that 



poo 

/ 1 + x 4 -x 2 {2 + x A ) l ' 2 dx 
Jo 



2 3 / 4 VHX3/4) 



5T(5/4) 

we see that the local energy to leading order is — Iqv(v / 1 3 ) 1 / 4 . For the one- 
component gas we should set v = pi 3 and for the two-component gas we 
should set v = (j) 2 £ 3 (see (|9.5[l ). After replacing the sum over boxes by an 
integral and at the same time replace the discrete Laplacian by a continuum 
Laplacian, as described above, we arrive at asymptotic lower bounds as in 
Theorems O and EP1 

There is one issue that we have not discussed at all and which played 
an important role in the treatment of the dilute gas. How do we know 



78 



E.H. LIEB, R. SEIRINGER, J. P. SOLOVEJ, AND J. YNGVASON 



the number of particles in each of the small cubes? For the dilute gas a 
super additivity argument was used to show that there was an equipartition 
of particles among the smaller boxes. Such an argument cannot be used for 
the charged gas. For the one-component gas one simply minimizes the energy 
over all possible particle numbers in each little box. It turns out that charge 
neutrality is essentially required for the energy to be minimized. Since the 
background charge in each box is fixed this fixes the particle number. 

For the two-component there is a-priori nothing that fixes the particle 
number in each box. More precisely, if we ignored the kinetic energy between 
the small boxes it would be energetically favorable to put all particles in one 
small box. It is the kinetic energy between boxes, i.e., the discrete Laplacian 
term in Theorem 19.51 that prevents this from happening. Thus we could in 
principle again minimize over all particle numbers and hope to prove the 
correct particle number dependence (i.e., Foldy's law) in each small box. 
This is essentially what is done except that boxes with very many or very 
few particles must be treated somewhat differently from the "good" boxes. 
In the "bad" boxes we do not prove Foldy's law, but only weaker estimates 
that are adequate for the argument. 



9.5. The Rigorous Upper Bounds. 

9.5.1. The upper bound for the two-component gas. To prove an upper 

(2) 

bound on the energy E y Q > (N) of the form given in Dyson's formula The- 
orem 19.21 we shall construct a trial function from the prescription in the 
Bogolubov approximation. We shall use as an input a minimizer <3> for the 
variational problem on the right side of ([9.4)1 . That minimizers exist can 
be easily seen using spherical decreasing rearrangements. It is however not 
important that a minimizer exists. An approximate minimizer would also 
do for the argument given here. Define <fr (x) = N 3 / 10 <^(N 1 / 5 x). Then again 
f <Pq = 1. In terms of the unsealed function (fr in ([9.5)1 . <fro(x) = A r_1 0(x). 

Let (f) a , a = 1, ... be an orthonormal family of real functions all orthogonal 
to (fro. We choose these functions below. 

We follow Dyson |D2j and choose a trial function which does not have a 
specified particle number, i.e., a state in the bosonic Fock space. 

As our trial many-body wave function we now choose 

^ = exp (-Aq + A ao+ + A aQ_) 

xHtl-^^expf- £ ^^ ee '<e<e')l°>> (9.9) 

a^O e.e'=±l a^O 

where a* a e is the creation of a particle of charge e = ±1 in the state cfr a , 
|0) is the vacuum state, and the coefficients Ao, Ai, . . . will be chosen below 
satisfying < A Q < 1 for a ^ 0. 

It is straightforward to check that ^ is a normalized function. 
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Dyson used a very similar trial state in D2 , but in his case the exponent 
was a purely quadratic expression in creation operators, whereas the one 
used here is only quadratic in the creation operators a* e , with q/0 and 
linear in ag_j_. As a consequence our state will be more sharply localized 
around the mean of the particle number. 

In fact, the above trial state is precisely what is suggested by the Bogol- 
ubov approximation. To see this note that one has 

(o 0± - A )^ = 0, and (a* + - a*_ + A a (a a+ - a a _)) * = 

for all Thus the creation operators for the condensed states can be 

replaced by their expectation values and an adequate quadratic expression 
in the non-condensed creation and annihilation operators is minimized. 
Consider now the operator 

00 A 2 

^ = ErV ww (9 - 10) 

a=l 

A straightforward calculation of the energy expectation in the state gives 
that 

(%f^HP*\ = 2A2 / u|(V0 o ) 2 + Tr(- /U A 7 ) 



N=0 



+2A 2 ,Tr ( K 



(k (7-^7(7 + 1))), 

where K, is the operator with integral kernel 

JC{x,y) = MVW-VI^MV)' (9- 11 ) 
Moreover, the expected particle number in the state ^ is 2A 2 , + 1^(7). In 
order for \E' to be well defined by the formula (|9,9|) we must require this 
expectation to be finite. 

Instead of making explicit choices for the individual functions 4> a and the 
coefficients A a , a / we may equivalently choose the operator 7. In defining 
7 we use the method of coherent states. Let X be a non-negative real and 
smooth function supported in the unit ball in M 3 , with J X 2 = 1- Let as 
before N~ 2/5 < £ <C N~ 1/5 and define Xi(x) = £~ 3/2 x(x/£). We choose 

JR 3 xIR 3 

where is the projection orthogonal to <po, 

0u,p( x ) = exp(ip • x)Xe(x-u), 

and 

//-mx 1 ( f + 167rA^/i-Vo(u) 2 A 

/KIpI) = o 



We note that 7 is a positive trace class operator, 7^0 = 0, and that all 
eigenfunctions of 7 may be chosen real. These are precisely the requirements 
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needed in order for 7 to define the orthonormal family cf) a and the coefficients 
A a for 

We use the following version of the Berezin-Lieb inequality |Be| IL3) . As- 
sume that £(t) is an operator concave function of M + U {0} with £(0) > 0. 
Then if Y is a positive semi-definite operator we have 

Tr (Y£( 7 )) > (2tt)- 3 j £(/(«, \p\)) V^YV^) dudp. (9.12) 

We use this for the function £(t) = yt{t + 1). Of course, if £ is the identity 
function then (|9.12|) is an identity. If Y = I then ()9.12j) holds for all concave 
i with £(0) > 0. 

Proving an upper bound on the energy expectation H9.11|) is thus reduced 
to the calculations of explicit integrals. After estimating these integrals one 
arrives at the leading contribution (for large Ao) 

2A§/i J (V0 O ) 2 + J f Uf + 2A^o(n) 2 ^) f(u, \p\) 

-_2A§0o(u) V/K \pW¥, \P\) + 1) dpdu 

= 2A^|(V0 o ) 2 -/o/(2Ag) 5 /Vf , 

where Iq is as in Theorem 19.21 

If we choose Ao = y N/2 we get after a simple rescaling that the energy 
above is iV 7 / 5 times the right side of (|9.4j) (recall that $ was chosen as the 
minimizer). We also note that the expected number of particles is 

2A 2 1 + Tr( 7 ) =iV + 0(iV 3 / 5 ), 

as iV — > 00. 

The only remaining problem is to show how a similar energy could be 

achieved with a wave function with a fixed number of particles N, i.e., how 

(2) 

to show that we really have an upper bound on Eq (N). We indicate this 
fairly simple argument here. 

We construct a trial function ^f' as above, but with an expected particle 
number N' chosen appropriately close to, but slightly smaller than N. More 
precisely, N' will be smaller than N by an appropriate lower order correc- 
tion. It is easy to see then that the mean deviation of the particle number 
distribution in the state is lower order than N. In fact, it is of order 
x/iV 7 ~ y/~N. Using that we have a good lower bound on the energy (n) 
for all n and that VP' is sharply localized around its mean particle number, 
we may, without changing the energy expectation significantly, replace 
by a normalized wave function ^ that only has particle numbers less than 

N. Since the function n 1— ► E^ (re) is a decreasing function we see that the 

(2) 

energy expectation in the state ^ is, in fact, an upper bound to E Q (N). 
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9.5.2. The upper bound for the one- component gas. The upper bound for the 
one-component gas is proved in a very similar way as for the two-component 
gas. We shall simply indicate the main differences here. We will again 
choose a trial state without a fixed particle number, i.e., a grand canonical 
trial state. Since we know that the one-component gas has a thermodynamic 
limit and that there is equivalence of ensembles |LN| , it makes no difference 
whether we choose a canonical or grand-canonical trial state. 

For the state (j)Q we now choose a normalized function with compact sup- 
port in A, that is constant on the set {x G A | dist(x,<9A) > r}. We shall 
choose r > to go to zero as L —* oo. Let us also choose the constant n 
such that ncpQ = p on the set where </>o is constant. Then n ~ pL 3 . 

Let again (p a , a = 1, . . . be an orthonormal family of real functions or- 
thogonal to As our trial state we choose, this time, 

* = [](! " A a) 1/4ex p(" A o/2 + Ao«S - E y«) 1°) ' ( 9 - 13 ) 

where a* is the creation of a particle in the state 4> a . We will choose \& 
implicitly by choosing the operator 7 defined as in (|9.1U|) . 
This time we obtain 




F - y\ 

+\ jj (p ~ Pi^) ~ x oM£) 2 ) \x - y\~ l (p - P y (y) - Ao<Ao(y) 2 ) dxdy 

AxA 

+Tr (- M A 7 ) + AgTr (k f 7 - VtCt + 1)) ) , (9-14) 



where P 7 (af) = j(x, x) and tC is again given as in (|9,11|) , We must show that 
we can make choices such that the first four terms on the right side above 
are lower order than the energy, and can therefore be neglected. 
We choose 

7 = 7£ = (2vr)- 3 1 fmKmmKM 

where e > is a parameter which we will let tend to at the end of the 
calculation. Here V^ as before is the projection orthogonal to 4>q and this 
time 

Tm) 2^(^ + 1^-^)1/2 J 

and 

9p(x) = \j np~ x exp(ip- x)4>q(x). 
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Note that np 1 <Pq(x) 2 is 1 on most of A. We then again have the Berezin- 
Lieb inequality as before. We also find that 

Pl {x) = (2vr)- 3 / f(\p\)dpnp- 1 Mx) 2 ( 1 + ( £ ~ lp ~ 1/4L ~ 1 )) 

= Aeip/fif^np^Mx) 2 (l + 0(e- l p~ l ' A L- 1 )) , 

where A £ is an explicit function of e. We now choose Ao such that A 2 , = 
n(l - Aep- 1 / 4 ^- 3 / 4 ), i.e., such that 

Xl<Pl{x) + Pl {x) = n0 o (f) 2 (l + 0(e- V 172 ^ 1 )) « P. 

It is easy to see that the first term in (|9,14jl is of order /5L 3 (rL) _1 and the 
fourth term in (|9.14|) is of order pL i {e~ 2 -\-pr 2 ). We may choose r, depending 
on L, in such a way that after dividing by pL 3 and letting L — > oo only the 
error e~ 2 remains. This allows choosing e -C p^ 1 ^ 8 . 

To estimate the second term in (|9.14(l we use Hardy's inequality to deduce 

^y^dxdy < 2(Tr 7 2 ) 1 / 2 Tr(-A 7 2 ) 1 /2 ) 
x-y\ 

and these terms can be easily estimated using the Berezin-Lieb inequality in 
the direction opposite from before, since we are interested now in an upper 
bound. The third term in ()9.14|) is controlled in exactly the same way as the 
second term. We are then left with the last two terms in (|9.14|) . They are 
treated in exactly the same way as for the two-component gas again using 
the Berezin-Lieb inequality. 
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